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INTRODUCTION 
Purpose 

One  ordinarily  defines  the  accidental  degeneracy  of  a  system  as 
that  which  does  not  follow  from  an  analysis  of  the  obvious  geometrical 
symmetry  of  the  system  in  configuration  space.   Frequently,  neverthe- 
less it  can  be  shown  that  there  is  a  higher  order  symmetry  group  in 
the  phase  space  of  the  system,  formed  with  the  help  of  certain  hidden 
symmetries,  which  is  adequate  to  account  for  this  accidental  degen- 
eracy.  A  well  known  example  of  a  system  possessing  this  property, 
and  one  of  the  first  to  be  investigated  historically,  both  classically 
and  quantum  mechanically,  is  the  non-relat i vist ic  Kepler  problem,  whose 
quantum  mechanical  analogue  is  the  problem  of  a  one-electron  atom. 

In  the  quantum  mechanical  case,  the  invariance  of  the  Hamiltonian 
with  respect  to  the  three  dimensional  rotation  group  accounts  for  the 
degeneracy  in  the  "m"  quantum  number,  but  not  for  that  in  the  "1" 
quantum  number. 

The  purpose  of  this  work  is  to  review  the  more  important  results 
obtained  on  the  accidental  degeneracy  of  harmonic  oscillators  and  the 
Kepler  problem  and  to  find  some  general  principles  with  which  the 
symmetry  groups  of  Hamiltonian  systems  can  be  found.   In  Chapter  1  a 
basic  mathematical  foundation  is  laid.   This  foundation  is  based  on 
the  fact  that  all  quadratic  functions  over  the  phase  space  of  a  system 
form  a  Lie  algebra.   On  this  basis  a  method  is  developed  which  enables 
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one  to  find,  rather  easily,  the  constants  of  the  motion  for  systems 
with  quadratic  Hami 1 tonians  or  for  those  which  can  be  made  quadratic 
through  a  canonical  transformation. 

This  method  is  then  applied  to  various  systems  in  the  succeeding 
two  chapters.   The  harmonic  oscillator  is  treated  in  Chapter  2,  both 
in  the  two  and  in  the  n  dimensional  case.   Chapter  3  deals  with  a 
general  Hamiltonian  which  may  be  interpreted  as  that  of  a  charged 
mass  point  in  a  plane  harmonic  oscillator  potential  and  uniform  mag- 
netic field.   Finally,  in  Chapter  4  the  Kepler  problem  is  discussed 
and  found  to  have  an  SU3  symmetry  besides  the  Rl^   generally  associated 
with  it. 

All  calculations  are  carried  out  classically.   The  results  ob- 
tained seem  to  be  valid  in  the  quantum  mechanical  case  except  as  noted. 
The  transition  to  the  quantum  mechanical  version  is  carried  out  in  the 
usual  manner  by  replacing  Poisson  bracket  relations  with  commutator 
braci<ets  and  replacing  functions  by  their  corresponding  quantum  mechani- 
cal operators. 

Historical  Sketch 

W.  Pauli  [1]  was  perhaps  the  first  to  associate  the  accidental  de- 
generacy of  the  Kepler  problem  with  the  less  well  known  second  vector 
constant  of  the  motion,  the  Runge  vector  R,  which  had  been  discussed 
earlier  by  W.  Lenz  [2].   Pauli  defines  this  vector  classically  as  lying 
in  the  plane  of  the  orbit  with  a  direction  from  the  force  center  to  the 


aphelion  and  a  length  equal  to  the  eccentricity  of  the  classical  orbit 
and  he  also  gives  the  quantum  mechanical  operator  definition.   The 
comutatlon  rules  between  the  components  of  this  vector  and  the  angular 
momentum  are  also  stated  explicitly.   As  a  consequence  of  the  the  com- 
mutation rules  he  shows  that  the  system  is  degenerate  in  "1"  also. 
However,  L.  Hulthen  [3]  observed  the  fact  that  these  commutation  rules 
are  the  same  as  those  of  the  generators  of  the  four  dimensional  rotation 
group.   0.  Klein  was  given  credit  for  first  recognizing  this  fact. 

V.  Fock  [k] ,     in  1935.  using  a  stereographic  projection  in  momentum 
space,  solved  the  Schrodinger  equation  in  momentum  space  and  showed 
that  the  solutions  are   spherical  functions  of  the  four  dimensional 
sphere  and  hence  that  the  symmetry  group  is  R^.   Later  V.  Bargmann  [5], 
in  a  discussion  of  Fock's  results,  showed  that  Fock's  group  was  gener- 
ated by  the  components  of  the  two  vector  constants  R  and  L,  the  angular 
momentum. 

It  was  about  this  time  that  0.  Laporte  [6]  obtained  a  differential 
operator  identical  to  the  Runge  vector  and  presented  a  method  for  ob- 
taining all  the  eigenfunct ions  of  the  hydrogen  atom  by  differentiation, 
similar  to  the  method  used  to  find  the  eigenfunct ions  of  L  by  using 
the  raising  and  lowering  operators  L.,.  and  L_.   He  made  extensive  use 
of  the  stereographic  parameters,  developed  by  him  and  G.Y.  Rainich  [7] 
about  the  same  time,  in  momentum  space.   The  wave  functions  were  found 
to  be  proportional  to  hyperspher leal  surface  harmonics. 

In  1939  and  1940  E.L.  Hill  and  his  student,  J.M.  Jauch  [8,  9,  10] 
published  two  articles  and  a  Ph.D.  thesis  dealing  with  accidental  de- 
generacy.  These  authors  considered  both  the  two  and  three  dimensional 


Kepler  problem  and  the  harmonic  oscillator  in  two  and  three  dimensions. 
The  Kepler  problem  was  found  to  have  R,  and  R^  symmetry  groups  in  two 
and  three  dimensions  respectively,  while  the  harmonic  oscillators  had 
the  symmetry  of  SU2  and  SUt,  the  groups  being  generated  by  the  constants 
of  the  motion  in  all  cases.   The  transformations  induced  by  the  gener- 
ators were  found  to  have  an  analogy  in  classical  mechanics,  namely  that 
they  represented  transformations  of  one  orbit  of  phase  space  into 
another  of  the  same  energy.   They  also  showed  the  existence  of  a  corre- 
spondence between  transformations  in  classical  mechanics  and  quantum 
mechanics. 

In  19^9,  A.W.  Saenz  [11]  wrote  a  Ph.D.  thesis  under  Laporte,  in 
which  the  symmetry  groups  explaining  the  accidental  degeneracy  of  vari- 
ous problems  were  found  by  reducing  the  problems  to  force  free  motion 
on  hyperspherical  surfaces  of  various  dimensions.  In  this  manner  he 
considered  the  Kepler  problem,  harmonic  oscillator  and  rigid  rotor. 

More  recently,  the  n  dimensional  isotropic  oscillator  has  been 
discussed  by  G.A.  Baker  [12].   He  found  the  synmetry  group  to  be  SU^ 
and  demonstrated  that  this  group  could  account  for  all  the  degeneracy 
found.   This  had  been  demonstrated  earlier  by  Y.N.  Demkov  [I3]  from  a 
slightly  different  point  of  view.   Whereas  Baker  attacked  the  problem 
by  finding  the  most  general  group  which  left  the  Hamiltonian  invariant, 
Demkov  used  the  generators  of  infinitesimal  transformations  which  com- 
muted with  the  Hamiltonian  and  then  showed  that  this  defined  an  SU^^ 
group.   In  a  later  article  Demkov  [l4]  demonstrated  the  equivalence  of 
his  and  Baker's  works. 


The  connection  between  accidental  degeneracy  and  hidden  symmetry 
was  confirmed  by  A. P.  Alliluev  [15]  in  1958.   He  considered  two  ex- 
amples, the  two  dimensional  oscillator  and  the  Kepler  problem  in  n 
dimensions,  which  reduced  to  Fock's  result  in  the  case  where  n  =  3. 
In  the  general  case,  the  hidden  symmetry  of  the  Kepler  problem  in  n 
dimensions  was  found  to  be  that  of  the  n  +  1  dimensional  hypersphere. 
The  one  dimensional  Kepler  problem  has  been  treated  by  R.  Loudon  [16] 
and  was  found  to  be  doubly  degenerate  in  all  levels  except  the  ground 
state,  an  apparent  violation  of  the  theorem  that  the  energy  levels  of 
one  dimensional  systems  are  nondegenerate.   Loudon  discusses  his  results 
in  the  light  of  this  theorem  and  shows  that  it  holds  only  when  the  po- 
tential is  free  of  singularities.   In  1959.  H.V.  Mcintosh  [17]  wrote  a 
qualitative  review  of  the  work  up  to  that  time  and  stated  the  results. 
The  most  recent  article  to  appear  is  one  by  Demkov  [18]  in  which  the 
concept  of  excessive  symmetry  groups  is  introduced.   He  gives  a  pre- 
scription for  finding  the  minimal  symmetry  group  and  as  an  example 
applies  it  to  the  anisotropic  two  dinnensional  oscillator  with  incom- 
mensurable frequencies. 

In  addition  to  the  works  cited  above  dealing  with  thestudy  of 
accidental  degeneracy  itself,  there  has  been  a  considerable  amount  of 
work  in  using  the  symmetry  groups  for  various  problems  and  applying 
them  to  the  many  body  problem,   particularly  to  the  shell  model  of 
the  nucleus.   One  of  the  first  of  these  was  by  J. P.  Elliott  [19,  20] 
who  undertook  to  study  collective  motion  in  the  shell  model  utilizing 
a  coupling  scheme  associated  with  the  degeneracy  of  the  three  dimensional 


harmonic  oscillator.   The  direct  result  was  a  classification  of  the 
states  according  to  the  group  SUo.   Following  this  treatment  was  a 
series  of  papers  by  V.  Bargmann  and  M.  Moshinsky  [21,  22]  on  the 
group  theory  of  harmonic  oscillators,  in  which  they  developed  a 
classification  scheme  for  states  of  n  particles  in  a  common  harmonic 
oscillator  potential.   They  were  primarily  concerned  with  finding  a 
scheme  which  would  exhibit  explicitly  the  collective  nature  of  the 
states.   Similar  work  was  also  done  by  S.  Goshen  and  H.J.  Lipkin 
[23,  m]   with  the  exception  that  the  particles  were  in  a  one  dimen- 
sional f>otentIal  well.   This  system  could  also  be  described  from  a 
collective  viewpoint  and  these  states  could  be  classified  using  the 
group  of  transformations. 

In  recent  times,  L.C  Biedenharn  has  extensively  studied  the 
Kepler  problem  and  its  symmetries,  both  in  the  non-relat i vist ic  and 
relativistic  cases.   In  the  non-relat i vist ic  [25]  case  he  has  obtained 
relations  for  the  representations  of  R.  exploiting  the  fact  that  Rk  is 
locally  isomorphic  to  Ro  x  R^.   The  non-relat ivist ic  problem  has  also 
been  considered  by  Moshinsky  [26],  where  use  has  been  made  of  the 
accidental  degeneracy  in  the  coulomb  problem  to  obtain  correlated  wave 
functions  for  a  system  of  n  particles  in  this  potential.   In  the  rela- 
tivistic case,  Biedenharn  [27]  and  Biedenharn  and  Swamy  [28]  have  ob- 
tained operators  which  are  analogous  to  the  angular  momentum  and  Runge 
vector  in  the  non-relat ivist ic  problem. 


CHAPTER  1 

MATHEMATICAL  FOUNDATIONS 

Introduction 

In  classical  and  quantum  mechanics  "quadratic"  Hamiltonians  play 
an  exceptional  role,  a  fact  which  is  traceable  in  the  end  to  the  fact 
that  the  Poisson  bracket  operation  satisfies  axioms  whose  expression 
is  tantamount  to  saying  that  the  Poisson  bracket  is  an  alternating  bi- 
linear functional  which  is  a  derivative  in  each  of  its  arguments. 

If  H^"'  designates  the  vector  space  of  homogeneous  polynomials 
of  degree  n,  whose  basis  may  be  chosen  to  be  composed  of  homogeneous 
monomials,  the  Poisson  bracket  operation  maps  the  cartesian  product 
^(m)  X  h(")  into  H^"^""^) ,  as  will  be  shown  later.   The  Poisson  bracket 
operation  is  a  bilinear  functional,  which  means  that  it  defines  a  linear 
transformation  when  either  argument  is  held  fixed,  and  so  it  always  has 
a  matrix  representation.   For  the  case  m  =  2,  one  has  a  mapping  from 
h(")  into  itself,  representable  by  a  square  matrix.   In  this  latter  case, 
the  extensive  lore  of  matrix  theory  is  available  to  help  discuss  the 
transformations.   In  particular,  the  mapping  may  be  described  in  terms 
of  its  eigenvalues  and  eigenvectors. 

Thus  the  mapping 

Tq(p)  =  (q.  p}.  (1.1) 

for  a  fixed  q  €  h'  ' ,    and  argument  p  e  h'"^,  can  be  expected  to  have 
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specialized  properties  which  become  all  the  more  significant  when 
n  =  2.   In  this  case  the  elements  of  H^  '  form  the  basis  of  a  Lie 
algebra  under  the  Poisson  bracket  operation.   In  any  event,  the  gen- 
eral mapping  T^  of  Eq.  (I.I)  yields  representations  of  this  algebra 
as  it  acts  on  polynomials,  homogeneous  of  various  degrees  as  shall 
be  shown  later  for  the  case  n  =  2.   Thus  the  allusion  to  the  excep- 
tional role  of  quadratic  Hamiltonians  refers  to  the  expectation  that 
the  algebraic  structure  of  the  Lie  algebra  will  have  its  reflection 
in  the  dynamical  properties  of  the  system  which  it  represents. 


Lie  Algebra  Defined  by  Poisson  Brackets 

The  set  of  2n  linearly  independent  coordinates  and  momenta  of  an 
arbitrary  physical  system  is  said  to  form  its  phase  space  0.   These 
variables  may  be  paired  in  such  a  way  that  they  can  be  indexed  by  fj, 
i  =  n,  ...I,  -1,  ...-n,  where  the  positive  indices  indicate  coordinates 
and  the  negative  ones  the  corresponding  conjugate  momenta.   The  Poisson 
bracket  of  any  two  functions  of  these  variables  is  defined  by 


''•"''Zl^srr-s^lTr-  <'-^' 


i  =  l 

It  follows  from  Eq.  (1.2)  that  the  Poisson  bracket  operation  satisfies 
the  axioms  of  a  Lie  algebra  [29,  3?]: 

(f.g)  =  -  [g.f)  (alternating  rule),      (1-3) 


{f,  ag  +  Bh}  =  a[f ,   g}  +  B(f,  h}  (linear),      {] .k) 

if.    ig.    h))  +  {g.  (h,  f})  +  {h,  {f,  g})  =  0  (Jacobi  identity). 

(1-5) 
The  Poisson  bracket  also  acts  like  a  derivative  in  that  it  satisfies 
the  relation 

{f.  gh)  =  g(f,  h}  +  [f,  g)h  .  (1.6) 

Although  all  second  degree  polynomials  form  a  Lie  algebra,  the 
linear  monomial  functions  of  the  coordinates  in  phase  space  itself 
also  form  a  linear  vector  space  0,  where  the  basis  may  be  chosen  to 
be  the  2n  fj.   Denoting  the  Poisson  bracket  operation  {  ,  }  as  the 
symplectic  inner  product  in  i> ,     it  is  noted  that  this  product  does 
not  satisfy  the  usual  axioms  for  an  inner  product  since  it  is  anti- 
symmetric under  interchange  of  the  arguments.   This  difficulty  may  be 
remedied  in  the  following  manner.   Using  the  definition  of  the  Poisson 
bracket  it  is  found  that 

1    i  =  -j  >  0 
(f..  f_.)  =0    i  ^  -j 

-I    i  =  -j  <  0  .  (1.7) 

From  this  it  is  observed  that  if  i  >  0,  the  dual  vector  f.  to  f.  may 
be  taken  to  be  f_j.  Similarly  for  i  <  0  the  dual  vector  is  taken  to 
be  -f  ..   Thus  the  dual  basis  to 

r=  (f^.  •••  f,.  f.,.  ...  f.n)  (l-8a) 
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is  then 

r*=  (-f.„.  .-.  -f...  fp  ...  V-    (1.8b) 


with  these  definitions  the  inner  product  may  be  defined  as 

(f.  g)  =  {f.  g}.  (1.9) 

where  f  and  g  are  linear  combinations  of  the  fj  and 


=  )^atfi  .  (1. 10) 


f 

i 


It  is  evident  then  that 

(f,.  fj)  =  (f,.  fj}  =6,j  (1.11) 


and 


(f.  f)  =  (f.  f)  =)^i^il^  •  <^-'2) 


The  space  <I>  also  has  the  peculiar  property  of  being  self-dual  in  a 
special  way,  namely  the  mapping  Tj:  is  a  linear  mapping  and  is  in  a  one 
to  one  correspondence  to  the  vectors  in  the  dual  space  $■*". 

Considering  tensor  products  of  $,  it  is  observed  that  (f,  g)  also 
defines  a  mapping  on  these  products.   If  one  defines 

h(")  =  <l>®<l><g) ®'i>  =  (^^,  (1-13) 

and  if  f  €  h("')  and  g  e  h(")  .  then  (f,  g)  6  h("^""^).   In  particular, 
when  m=  2  i.e.,  f  e  Q  s  h'  ' ,  a  mapping  is  obtained  which  is  represent- 
able  by  a  square  matrix. 


n 


Such  a  mapping  may  be  expected  to  have  eigenfunct ions  and  eigen- 
values.  Let  h  be  a  fixed  element  of  Q  (a  quadratic  Hamiltonian)  and 
let  g;  and  \-    be  the  eigenfunct ions  and  eigenvalues  such  that 

(h.  g;)  =  Th(g;)  =  ^jg;.  (1.14) 

where  the  gj  are  taken  to  be  linear  I.e.  g;  €  H^  '.   If  h  is  represent- 
able  as  a  normal  matrix,  the  gj  form  a  basis  for  the  space  h'  ',  and 
are  simply  linear  combinations  of  the  f|.   In  particular,  since  the 
gj  €  <&  s  h'  ',    they  form  a  basis  for  all  the  tensor  spaces  as  well 
where  all  possible  products  of  the  gj  are  taken  i.e., 


((gi))  =  H<') 

((9,9^))  =  H^2) 


((gj....  9j)    =  H<").  (1.15) 


On  account  of  the  derivative  rule  Eq.  (1.6)  the  following  relation  is 
satisfied 

{h,  gjgj)  =  g;(h,  g j )  +  {h,  gj)  gj 

=  (K;  +  Xj)  gjgj,  (1.16) 

and  hence  the  monomials  in  gj  are  eigenfunct ions  of  h  in  every  H^"' 
and  form  a  basis  even  in  the  case  of  degeneracy,  for  all  possible  eigen- 
funct ions. 

In  cases  where  the  operator  Tf  is  normal,  a  complete  set  of  or- 
thogonal eigenvectors  exists  which  will  be  denoted  g..   Hence,  if  the 
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eigenfunctions  are  indexed  such  that  gj  =  g_j  then  the  eigenvectors  of 
any  normal  operator  satisfy  the  rule 

(gj.  9j}  =  &j,.j.  (1.17) 

The  eigenvalues  \|  corresponding  to  the  eigenvectors  gj  of  the  normal 
operator  h  occur  in  negative  pairs.   This  follows  from  a  consideration 
of  the  Jacobi  identity 

lh,{g_.,  g.)}  +  {g_j,  {g.,  h})  ^■  {g..{h.  g_.))  =  0. 

thus 

(h,  1}  +  (-\.)  [g_..  g.}  +  \_.  {g.,  g_.)  =  0. 

and 

(^-i  +  ^0   {gj.  g-i)  =  0,  (1.18) 

and  hence,  \.   =   -\    .    since  {q.,  q  .}  /  0. 


Since  h  e  H^  ' ,    h  may  be  written  as 


n 

=  )  C. .g.g. , 


h  =  )  C,.g.g..  (1.19) 


i,j=-n 


•  =  J 


where  the  g;  are  eigenfunctions  of  h  belonging  to  eigenvalues  \j  and 

are  indexed  so  that 

K^   =   -\_i  .  (1.20) 

Assuming  that  the  eigenvalues  of  h  are   distinct  and  using  the  identity 
{h,  h)  =  0  one  has,  with  the  aid  of  Eq.  (1-19) 
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0  =\   C..[h.  g.g.} 


«.J 
i2y 


\  C..(\.  +  \.)  g.g..         (1.21) 


'.J 


Hence  each  coefficient  must  vanish  and  either  Cj  j  =  0  or  {\-    +  >.:)  =  0. 
However  X.?  =  -\j  only  if  j  =  -i.   Therefore 

n 

h=^  Cj^.j  g;g.;.  (1.22) 

i=i 


To  further  evaluate  these  coefficients  recall  that 

{h,  g^)  =  \^g^.  (1.23) 

thus  A 

H9k  =2.  ^!.-'  ^^'^-''   ^'<^ 

i«l 
n 

Using  the  fact  that  the  coefficients  C..  are  symmetric  in  the  indices 
the  coefficient  of  g^  on  the  right  side  is 


\=Ck,-k^^-k'  ^k^'  <'-25) 


]k 


and  hence 

C.  .^^ ('-26) 


The  final  conclusion  is  that 


h  =  ■ 


-Iji^'-''-'-  <'■"' 


l=i 


I    I 


Consider  now  the  case  where  h  may  have  degenerate  eigenvalues. 
Using  the  general  form  for  h,  Eq.  (1.19),  and  calculating  {h,  h)  as 
before  Cj j  =  0  unless  X.  =  -\. .   However  in  the  present  instance  it  is 

possible  to  have  X.  =  \  ,  k  ?^  m  and  hence  Xi  =  -X  „.   Hence 

"^  km  k  -m 


g.g      .  (1.28) 

, -m   ^k^-m  * 


n 

n 

i  =  l 

k,m=l 

X  =X 
k     m 

k  ?^  m 

Calculating  once  again    {h,    g.)    results    in   the   following   equations; 

n  n 

^^j  =_X    S.-i    ^^i^-i-  9.)  *^    C,  _^  (g,g_^.  g.)         (1.29) 

i=l  k,m=l 

X,  =X 
k     m 

k  ?^  m 


II 

=  ^    C...    [g,    {g...   g.)  +  g.j    (gj.   g.)] 
i=l 
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+  )  C,    [g, {g  ,  g.)  +  g  [g. ,  g.}]  . 
/  ^   k,-m  ^k  ^-m  ^j    ^-m^k*  ^j 

k,nT=l 

^k=Nn 
k  it  m 

Consider  the  first  term  in  the  second  sum.   If  m=j  then  k?^  implies 
k.^}    and  hence  C.  _.  must  be  zero.   Therefore,  as  before 


h  = 
i  =  I 


II 

y  -r-^ rgig-i-  (1-30) 


It  should  also  be  noted  that  the  transformation  to  eigenvector 
coordinates  for  h  is  a  canonical  transformation  since  the  eigenvectors 
satisfy  (gj,  gj)  =  Sj  _.,  which  is  a  necessary  and  sufficient  condition 
for  g.  and  g  .  to  be  canonical  variables. 

Consider  now  the  canonical  form  of  a  non-normal  quadratic  operator 
h,  of  the  Jordan  canonical  form.   Again 


n    n 
"=11   C,.g,g..  (,.3.) 


i  =  -n  j«=-n 


where  the  gj  are  the  principal  vectors  of  h  and  indexed  such  that  the 
same  orthogonality  conditions  hold  as  before.   Using  the  property  that 


(h,  g.)   =  g;^,.  (1.32) 


it  is  found  that 
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i=-n  j=-l 
i  i  J 

=  )   )  C   [g  {g.,  g  }  +  g.{g.,  g,  )] 

l=-n  j=-n 
i  i  J 

'Ck+l,-k3k+l     '<>0 

.C-k.k+l  Sk+l     •<"  0  .  (1.33) 


Hence  all  C: ;  =  5.  _.  ,  and 

n 

h=)^  gig.i^,  .  (1.3^) 

i=l 

The  sympletic  norm  also  produces  certain  symmetries  and  anti- 
symmetries.  For  example,  the  dual  mapping  R(fj),  which  maps  the 
coordinates  into  momenta  and  the  momenta  into  the  negatives  of  the 
coordinates , 

i  >  0  (1.35) 

R(f  .)  =  -f. 

preserves  the  inner  product;  i.e. 

(R(fi).  R(fj)}  =  {f;.  fj)  .  (1-36) 

On  the  other  hand  the  mapping  P(f:)  which  simply  exchanges  coordinates 
and  momenta 
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P(fi)  =  f-i.  (1.37) 

or  the  time  reversal  mapping  T(fj) 

Kfj)  =  f, 

i  >  0  (1-38) 

T(f.i)  =  -f.i 

reverses    the  sign   of   the  symplectic  norm 

{P(fj).    P(fj))   =   -    {f,.    fj] 

(T(f.).   T(f.))   =   -    (f..    f.}    .  (1.39) 

These  operators  may  be  expressed  in  terms  of  the  eigenvectors,  g.  of 

the  normal  operator  h,  previously  defined.   Recall  also  that  the  gr 

form  a  set  of  canonical  coordinates  and  momenta.   Using  the  convention 

that  positive  indices  represent  coordinates  and  negative  indices  the 

corresponding  canonical  momenta,  these  operators  have  the  following 

canonical  forms: 

n 


i=i 

2(g. 

,9,) 

n 

-I 

i  =  l 

g;g;  ■ 

•  9-i9.i 

2(9i 

9.,) 

n 

-I 

i=l 

919. 

-i 

(l.ill) 


{\M1) 
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These  operators  have  the  following  effects  on  eigenvectors  of  h: 

r^k       •<  >  0 
1^(9^)  "  (T.  g^)  =  j  (1.43) 

Tp(g,^)  =  {P.g,^}  =  g_,^        all  k,        (1.^) 

f  g.k         •<  >  0 

L-9-k       ^ "  °  • 

The  composite  of  any  two  of  the  operations  on  an  eigenfunct ion  is 
equivalent  to  operating  with  the  third  operator: 


M-k   ^^  0] 


VV9,)).9.,  =  yg,)  (I.V) 


It  is  also  interesting  to  note  that  the  three  operators  satisfy  the 
following  Poisson  bracket  relations  among  themselves: 

(P.  T)  =  2R  (1.^9) 

(R.  T)  =  2P  (1.50) 

{p.  R)  =  2T  (1.51) 

The  composite  of  the  normal  operator  h  with  each  of  these  operators 
on  the  eigenvectors  is  also  worth  noting.   It  is  relatively  easy  to  show 
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that 

and  hence  the  matrices  representing  h  and  T  commute.   However  one  finds 
that  P  and  R  anticommute  with  h,  i.e. 


VW--VVV>  ('•") 


and 


T  (T  (g  ))  =  -  T  (T  (g  )).  (1.54) 

R  h  k        h  R  k 


The  existence  of  either  one  of  these  last  equations  is  sufficient  to 
insure  that  the  eigenvalues  appear  in  negative  pairs  [30]. 

In  the  next  two  chapters  the  method  described  above  will  be  applied 
to  harmonic  oscillator  Hamiltonians  in  two  and  n  dimensions  and  to  vari- 
ations on  the  oscillator  Hamiltonian. 


CHAPTER  2 

THE  HARMONIC  OSCILLATOR 

Introduct  ion 


The  accidental  degeneracy  and  symmetries  of  the  harmonic  oscillators 
has  been  discussed  by  several  authors.   A  particularly  thorough  account 
has  been  given  in  the  thesis  by  J.M.  Jauch  [9],  in  which  he  considers  the 
isotropic  oscillator  in  two,  three  and  n  dimensions,  and  finds  the  sym- 
metry group  to  be  SU^,  i.e.,  the  Hamiltonian  is  invariant  under  transfor- 
mations induced  by  the  generators  of  the  group.   A  short  discussion  of 
the  anisotropic  oscillator  in  two  dimensions  with  commensurable  frequen- 
cies is  also  given.   The  main  results  of  this  work  are  also  contained  in 
the  notes  of  E.L.  Hill  [33]-   More  recent  discussions  of  the  n  dimensional 
isotropic  oscillator  have  been  given  by  G.A.  Baker  [12]  and  Y.N.  Demkov 
[13,  14.  18]. 

The  plane  Isotropic  Oscillator 

In  units  such  that  the  mass  and  spring  constant  are  unity,  the 
Hamiltonian  of  the  plane  isotropic  harmonic  oscillator  is 

H  =  (P?  +  Pj  +  x^  +  y2)/2  .  (2.1) 

X     y 

With  respect  to  the  basis  (x,y,P  ,P  ),  the  matrix  representation  of  H 
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2] 


as  an  operator  under  the  Poisson  bracket  relation  is. 


0 

0 

1 

0 

0 

0 

0 

1 

I 

0 

0 

0 

0 

-1 

0 

0 

(2.2) 


This  representation  is  found  by  calculating  the  effect  of  H  on  each 
member  of  the  basis  under  Poisson  bracket.   With  the  convention  that 
X  is  a  column  matrix  with  a  I  in  the  first  row,  y  with  a  1  in  the 
second  row  and  so  on,  the  matrix  representation  of  H  is  then  easily 
constructed.   The  eigenvalues  of  the  matrix  H  are 


\  =  ±  i, 


(2.3) 


each  root  appearing  twice.   The  normalized  eigenvectors  associated 
wi  th  \  =  +  i  are 

a  =  (1//2)(P^  -  ix)  =  1//2  f  o|    b  =  (I//2)(P^  -  iy)  =  lATa  ^ 

0. 


(2.M 


and   those  associated  with  X  =   -    i    are 


a*  =   (lAr2)(P^  +    ix)   =    ]/J'2   (O 

\0l 


b*=    (lAr2)(Py   +    iy)    = 


iAr2 


(2.5) 
where  the  column  matrices  are  the  matrix  representations  of  the 
functions  on  the  left.   One  may  easily  verify  this  representation, 
for  example 
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{H.  a]  =  (H.  P^  -  ix}  =  i(P^  -  ix)  (2.6) 

while  multiplying  the  matrix  in  Eq.  (2.2)  into  the  column  matrix  repre- 
senting a  also  gives  i  times  the  same  column  matrix.   The  quantum  me- 
chanical analogues  of  these  four  quantities  are  the  well  known  raising 
and  lowering  operators  for  the  harmonic  oscillator.   With  these  four 
quantities  ten  linearly  independent  quadratic  monomials  can  be  con- 
structed which  are  eigenfunct ions  of  H  and  have  as  eigenvalues  0,  +  2i 
or  -  2i.   Only  four  of  these  are  of  interest  here,  namely  those  with 
eigenvalue  zero.   They  are  the  four  linearly  independent  constants  of 
the  motion  (since  they  commute  with  H  under  Poisson  bracket)  listed 
be  1 ow : 

aa*   =  (p2  +  x2)/2  (2.7) 

bb*=  (p2  +  y2)/2  (2.8) 

ab"'-  =  [(P^Py  +  xy)  +  i(yP^  -  xP^)  ]/2  (2.9) 

a*b  =  [(P^Py  +  xy)  -  i(yP^  -  xPy) ]/2  .         (2.10) 

For  the  purpose  of  physical  interpretation  it  is  more  convenient  to 
deal  with  the  real  and  imaginary  parts  of  these  constants  separately, 
as  well  as  to  separate  the  Hamiltonian  from  them.   Accordingly  the 
following  four  quantities  are  introduced: 

H  ■=  aa'^   +  bb*  =  (P^  +  p2  +  x^  +  y2)/2         (2.11) 
X    y 

D  =  aa'-'^  -  bb'-^  =  (p2  +  x^  -  pj  -  y2)/2         (2.12) 
X  y 
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L  =  i(a*b  -  ab*)  =  yP^  -  xP, 


K  =  a""b  +  ab'^  =  xy  +  P^P 


y  ■ 


(2.13) 
(2.14) 


These  four  constants  have  the  following  simple  interpretations. 

a)  H,  the  Hamiltonian,  is  the  total  energy  of  the  system. 

b)  D  is  the  energy  difference  between  the  two  coordinates. 

c)  L  is  the  angular  momentum  of  the  system  and  the  generator  of 
rotations  in  the  xy  plane. 

d)  K  is  known  as  the  correlation  and  is  a  peculiar  feature  of  the 
harmonic  oscillator.   As  a  generator  of  an  infinitesimal  con- 
tact transformation,  it  generates  an  infinitesimal  change  in 
the  eccentricity  of  the  orbital  ellipse  while  preserving  the 
orientation  of  the  semiaxes,  and  preserving  the  sum  of  the 
squares  of  their  lengths.   The  energy  of  a  harmonic  oscillator 
depends  only  on  the  sum  of  the  squares  of  the  semiaxes  of  its 
orbital  ellipse,  which  remains  constant  under  such  a  transfor- 
mat  ion. 

The  set  of  functions  {K,  L,  D}  is  closed  under  the  Poisson  bracket 
operation.   Explicitly,  their  Poisson  bracket  table  is 


(2.15) 


K 

L 

0 

K 

0 

2D 

-2L 

L 

-2D 

0 

2K 

0 

2L 

-2K 

0 

Aside  from  the  factor  2  these  are  the  Poisson  bracket  relations  of  the 


2k 


generators  of  the  three  dimensional  rotation  group,  or  the  three  com- 
ponents of  the  angular  momentum  in  three  dimensions.   Since  the  con- 
figuration space  of  the  two  dimensional  isotropic  oscillator  is  appar- 
ently only  two  dimensional  and  has  only  rotations  in  the  xy  plane  as 
an  obvious  symmetry,  the  occurrence  of  the  three  dimensional  rotation 
group  is  rather  anomalous. 

By  performing  a  series  of  geometrical  transformations  on  a,  a*,  b 
and  b'  one  can  explicitly  demonstrate  the  spherical  symmetry  of  the 
system.   This  series  of  transformations  is  sometimes  called  the  Hopf 
mapping. 

The  Hopf  Mapping  [31] 

Perhaps  the  best  analytic  representation  of  the  Hopf  mapping  is 
obtained  by  introducing  polar  coordinates  (\,  t,  a,  p)  for  the  complex 
variables  a  and  b.   Explicitly,  these  are 

a  =  (\e-'P  cos  t)//2  (2.16) 

b  =  (\e-'<'  sin  i)Ar2,  (2.17) 

where  \  does  not  denote  the  eigenvalue  of  H  as  previously.   The  first 
step  in  the  Hopf  mapping  is  to  form  the  ratio  to  defined  as  the  quotient 

u)  =  X^  a/b  .  (2.18) 

The  formation  of  this  ratio  may  be  regarded  as  a  gnomonic  projection  of 
the  four  dimensional  space,  regarded  as  having  two  complex  dimensions, 
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onto  a  two  dimensional  space  (Fig.  2.1),  having  one  complex  dimension. 
This  maps  the  point  (a,b)  lying  on  the  two  dimensional  complex  circle 
of  radius  \   into  the  point 

X^  a/b=  \2  e'^*^^^  cot  T  (2.19) 

on  the  complex  line. 

The  second  step  is  to  regard  the  complex  point  u  =  \2  gi (a-p)  ^q^  ^ 
as  not  lying  on  a  complex  line;  but  rather  as  a  point  in  a  real  two 
dimensional  plane.   One  now  performs  an  inverse  stereographic  projection 
onto  the  surface  of  a  three  dimensional  sphere  whose  south  pole  is  tan- 
gent to  the  plane  at  its  origin.   By  choosing  the  origin  of  the  three 
dimensional  space  to  be  at  the  center  of  the  sphere,  a  point  on  the 
surface  of  the  sphere  will  be  specified  by  giving  the  radius  r,  the 
azimuth  4>  and  the  colatitude  9.   For  convenience  the  sphere  is  chosen 
to  have  diameter  \   .      The  azimuth  is  measured  in  a  plane  parallel  to  the 
original  plane  and  hence 

4>  =  o  -  p  .  (2.20) 

Thus  to  determine  the  location  of  the  projected  point  P  on  the  sphere, 
only  the  colatitude  Q    is  left  to  be  determined.   From  Fig.  2.2  one  has 


tan  (f  -|)  =  cot  T  (2.21) 


but 


tan  (f  -|)  =  cot  I  (2.22) 

hence 

e  =  2t  .  (2.23) 
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Thus  the  coordinates  of  the  point  P  projected  onto  the  surface  of 
the  sphere  from  the  point  w  =  X^  e'(^'P)  cot  t  are 

e  =  2t  (2.24) 

<t>  =  a  -  p  (2.25) 

r  =  \2/2  .  (2.26) 

Expressing  the  quantities  H,  K,  L  and  D  in  Eq.  (2.4)  in  terms  of 
these  new  coordinates  gives 

H  =  \^n  (2.27) 

K  =  (X^/2)  sin  0  cos  <l>  =  X             (2.28) 

L  =  (\^/2)  sin  e  sin  *  =  Y             (2.29) 

D  =  (\2/2)  cos  e  =  Z  .  (2.30) 

Hence,  K,  L  and  D  simply  determine  a  point  on  the  surface  of  this  sphere. 
It  should  be  noted  that  K  generates  infinitesimal  rotations  about  X,  L 
about  Y,  and  D  about  Z. 

Since  these  three  constants  of  the  motion  are  then  directly  the 
coordinates  of  the  fixed  point  corresponding  to  the  orbit  of  the  oscil- 
lator, it  is  apparent  that  they  generate  rotations  of  the  sphere,  which 
transform  one  orbit  into  another  of  the  same  energy. 
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Canonical  Coordinates 

Whenever  two  functions  satisfy  the  Poisson  bracket  relation 

U.  B}  -  HB  (2.31) 

where  p.  i s  a  scalar,  a  corresponding  equation, 

(A,  (1/n)  In  B)  =  1  ,  (2.32) 

may  always  be  written,  which  is  the  rule  satisfied  by  a  coordinate  and 
momentum  which  form  a  canonical  ly  conjugate  pair.   The  first  of  these 
two  equations  simply  states  that  B  is  an  eigenfunct ion  of  A  with  respect 
to  the  Poisson  bracket  operation.   Consequently  it  is  to  be  expected 
that  such  eigenfunct ions  be  a  ready  source  of  canonical  variables.   A 
system  in  which  the  Hamiltonian  is  one  of  the  momenta  can  be  obtained 
in  the  present  case.   Consider  the  following  diagram  in  which  the  two- 
headed  arrows  imply  that  the  Poisson  bracket  of  the  two  functions  is 
zero  and  the  one-headed  arrows  imply  that  the  Poisson  bracket  is  +  1 
in  the  direction  indicated: 


H  =  a  a  +  b  b'^    < >  D  =  a  a"^  -  b  b 


(lAi)  In  {aWab)< >  (lA')  In  (ab"Va*b) . 

(2.33) 
If  these  four  quantities  are  rewritten  in  terms  of  the  four  parameters 
introduced  in  the  Hopf  mapping  and  p  +  o  is  denoted  by  ^,  then  Eq.  (2.33) 
becomes 
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^/2  <r 


>  D  =  (X  /2)  cos  9 


■^  4)/2 


(2.3^) 


Hence  the  Hamiltonian  and  the  energy  difference  are  the  two  momenta 
and  the  angles  Y/2  and  4>/2  are  the  corresponding  canonical  coordinates. 
Since  ^/2  is  conjugate  to  H,  "i    increases  linearly  with  time. 

The  final  result  of  this  transformation  is  to  map  an  entire  orbit 
into  a  single  point.   Since  the  radius  of  the  image  sphere  determines 
the  energy,  the  only  Information  lost  as  a  consequence  of  the  mapping  is 
the  phase  of  the  point  representing  the  harmonic  oscillator  in  the  great 
circle  comprising  its  orbit.   The  lost  phase  may  be  recovered  by  attach- 
ing a  flag  to  the  point  representing  the  orbit,  which  will  then  rotate 
with  constant  velocity.   The  missing  angle  Y  may  then  be  taken  as  the 
angle  which  the  flag  makes  with  its  local  meridian.   In  this  way  the 
motion  of  the  harmonic  oscillator  in  its  phase  space  is  shown  to  be 
equivalent  to  the  motion  of  a  rigid  spherical  rotator,  and  in  fact  the 
parameters  4),  0  and  ¥  are  just  the  Euler  angles  describing  the  orien- 
tation of  the  rotor.   These  same  parameters  may  also  be  used  for  the 
description  of  spinors,  and  a  particularly  lucid  geometrical  interpre- 
tation of  this  transformation  may  be  found  in  a  paper  of  Payne  [32] 
describing  various  methods  of  representation  of  two-component  spinors. 
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The  Anisotropic  Harmonic 
Oscillator  in  Two  Dimensions 


The  anisotropic  oscillator,  both  with  commensurable  frequencies 
and  with  non-commensurable  frequencies  has  been  discussed  by  Jauch  and 
Hill  [10]  and  also  by  Hill  [33]-   These  authors  were  able  to  show  that, 
in  the  case  of  commensurable  frequencies,  the  symmetry  group  of  the 
system  is  again  the  three  dimensional  rotation  group.   However,  they 
were  unable  to  come  to  any  conclusion  concerning  the  case  with  incom- 
mensurable frequencies.   Recently  R.L.  Hudson  has  been  investigating 
the  incommensurable  case  also. 

It  can  be  shown  though,  that  even  with  non-commensurable  frequen- 
cies the  symmetry  group  is  also  a  three  dimensional  rotation  group. 
The  proof  is  sketched  below  and  follows  closely  that  used  in  the  iso- 
tropic case,  except  where  noted. 

If  units  are  chosen  such  that  the  frequency  in  the  x  coordinate  is 
unity,  the  Hamiltonian  is 

H  =  (p2  +  p2  +  x2  +  u)2y2)/2  ,  (2.35) 

where  oj  is  the  frequency  in  the  y  coordinate  and  is  an  irrational  number. 
The  eigenfunct ions  and  eigenvalues  of  H  under  the  Poisson  bracket  opera- 
tion are  listed  below: 


Private  communication. 
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Eigenfunction 

Eigenvalue 

a  =  (lAr2)(Px  -  ix) 

a*  =  (1/^2)(P^  +  ix) 

-i 

b  =  (l//2)(Py  -  iwy) 

iu) 

b*  =  (lAr2)(Py  +  icjy) 

-iU) 

(2.36) 

Although  the  most  obvious  four  constants  of  the  motion  are  aa  ,  bb  , 
(a)*^*,  and  (a")'%,  a  more  convenient  set  is  the  following  collection  of 
independent  combinations  of  them: 


H  =  aa*  +  bb* 


D  =  aa   -  bb 


K=  [(a)V*  (a*)"b]/(aa*)(^-')/2 
L=  i[(a)V  -  (a*)'^]/(aa*)('^l)/2 


(2.37) 
(2.38) 
(2.39) 
(2.40) 


Again  the  set  {K,  L,  D}  is  closed  under  the  Poisson  bracket  operation 
and  forms  a  Lie  group,  the  three  dimensional  rotation  group.   The 
Poisson  bracket  table  is  listed  below: 

L         D 


K 

K 

0 

L 

-2a)D 

D 

2u)L 

2(00 
0 
-2(iJK 


-2u)L 
2ajK 
0 


(2.41) 


A  mapping  to  polar  coordinates  as  in  the  isotropic  case  can  once 
again  be  performed,  however  it  is  not  a  Hopf  mapping.   The  form  of  the 
final  equations  for  the  constants  is  the  same  as  in  Eqs.  (2.2?)  to 
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(2.30),  however,  for  the  anisotropic  case 


♦  »  Cijp  -  o  . 


(2.42) 


A  set  of  functions  which  proves  to  be  convenient  to  use  as  the  new 
canonical  coordinates  is  diagrammed  below  with  the  same  convention  as 
before: 


In  -* * < *•  In  ■* '— X 

kiui  (a)*^         4iu)    (a)'^* 


(2.i*3) 


If  the  parameters  from  the  mapping  are  substituted  into  Eq.  (2.43) 
this  table  becomes  identical  to  Eq.  (2.34)  with  the  conditions  that  4> 
Is  defined  by  Eq.  (2.42)  and 


^  =  up  +  a 


(2.44) 


It  is  rather  interesting  to  note  that  the  entire  effect  of  the 
possible  incommensurability  of  the  frequencies  in  the  two  coordinates 
is  completely  absorbed  by  the  power  mapping  (a)*^,  which  is  well  defined 
as  exp  (u)  In  a)  in  the  case  of  an  irrational  u).   In  this  latter  case 
one  has  a  mapping  with  infinitely  many  branches  rather  than  a  finite 
number  as  in  the  case  of  a  rational  w.   It  is  of  course  on  this  account 
that  the  orbit — a  Lissajous  figure  in  the  x-y  plane--would  then  be  open, 
and  would  not  correspond  to  periodic  motion. 
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Conmensurabil ity  and  Constants  of  the  Motion 

In  reviewing  the  analysis  of  the  theory  of  the  plane  harmonic 
oscillator,  it  is  probably  appropriate  to  draw  attention  to  a  number 
of  statements  made,  both  in  the  literature  of  classical  mechanics, 
and  the  lore  of  accidental  degeneracy.   First  of  all ,  the  existence  of 
accidental  degeneracy  is  often  considered  to  be  allied  with  the  exist- 
ence of  bounded,  closed  orbits.   When  this  argument  is  used,  it  is 
applied  to  the  existence  of  algebraic  integrals  of  the  motion,  and  in- 
deed It  should  be  noted  here  that  as  long  as  the  frequencies  of  the  two 
coordinates  are  commensurable,  bounded  closed  orbits  result  and  the  in- 
tegrals are  algebraic,  as  required  by  the  theory.   However,  the  above 
results  also  verify  something  which  has  been  generally  known  [3^].  namely 
that  non-algebraic  integrals  can  exist,  and  as  would  be  expected,  are 
associated  with  space-filling  motion. 

Another  point  regards  the  separability  of  the  Ham i 1  ton- Jacob i 
equation  in  several  coordinate  systems.   It  is  pointed  out  that  if  a 
Hamilton-Jacobi  equation  is  separable  in  more  than  one  coordinate  system, 
it  is  necessarily  degenerate  [35]-   The  argument  usually  given  concerns 
the  fact  that  orbits  with  incommensurable  frequencies  in  one  coordinate 
system  will  be  space-filling  curves,  the  boundaries  of  which  are  co- 
ordinate arcs.   Thus,  in  whatever  coordinate  system  the  motion  is  de- 
scribed, the  system  of  separation  is  uniquely  defined  as  that  one  form- 
ing an  envelope  of  the  various  orbits.   Of  course,  the  option  left  in 
the  enunciation  of  this  theory  is  the  fact  that  changes  of  scale  may 
still  be  made,  in  which  the  new  coordinates  are  functions  each  of  only 
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one  of  the  old  coordinates,  and  not  several  or  all.   So  it  is  in  the 
polar  mapping,  where  one  forms  a  ratio  of  powers  of  functions  of  only 
a  single  coordinate.   However,  this  change  of  "scale"  seems  to  eliminate 
the  problem  of  the  incommensurable  frequencies  at  the  same  time,  and 
thus  circumvents  the  "proof 'of  uniqueness.   This  result  has  only  been 
secured  at  a  price,  namely  the  inverse  of  the  transformation  is  infin- 
itely multiple  valued  [36].   Otherwise,  one  could  easily  reduce  the 
anisotropic  oscillator  to  a  static  point  on  the  Riemann  sphere,  return 
from  this  to  an  isotropic  oscillator,  and  separate  in  whatever  system 
desired. 

These  considerations  are  not  unrelated  to  the  problem  of  finding 
a  Lie  group  generated  by  the  constants  of  the  motion.   In  this  respect, 
a  Lie  group  has  always  been  regarded  as  being  more  important  quantum 
mechanically  than  classically,  because  its  irreducible  representations 
would  be  associated  with  accidental  degeneracies  of  the  Schroedinger 
equation.   Since  the  anisotropic  harmonic  oscillator  has  no  accidental 
degeneracies  in  the  case  of  incommensurable  frequencies,  it  is  some- 
what surprising  to  find  that  the  unitary  unimodular  group  is  still  its 
symmetry  group,  at  least  classically. 

In  the  case  of  commensurable  frequencies,  proportional  say  to  m 
and  n,  the  operator  a'"b'''"  commutes  with  the  Hamiltonian,  and  is  inter- 
preted as  creating  m  quanta  in  the  first  coordinate  and  annihilating  n 
quanta  in  the  second  coordinate.   Thus,  m  quanta  in  one  coordinate  are 
equal  to  n  quanta  in  the  second  coordinate,  the  energy  is  restored  to 
its  original  level,  and  the  system  is  left  in  another  eigenstate. 
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Thus,  it  is  possible  to  expect  and  explain  degeneracy  in  the  case  of 
commensurable  frequencies. 

In  the  classical  incommensurable  case,  the  symmetry  group  still 
exists.   After  all,  the  total  energy  of  the  oscillator  is  the  sum  of 
the  energies  in  each  coordinate,  which  is  reflected  in  the  fact  that 
it  is  the  sum  of  the  squares  of  the  semiaxes  of  the  orbit  (A  Lissajous 
figure,  like  an  ellipsoid,  has  principal  axes),  and  not  the  individual 
semiaxes,  which  determine  the  energy.   Thus,  the  dimensions  of  the 
bounding  rectangle  of  a  space-filling  Lissajous  figure  may  be  changed 
subject  to  this  constraint.   This,  in  addition  to  an  adjustment  of  the 
relative  phases  in  the  two  coordinates,  is  the  effect  of  the  constants 
of  the  motion  acting  as  the  generators  of  infinitesimal  canonical  trans- 
formations. 

The  problem  of  carrying  this  result  over  to  quantum  mechanics  seems 
to  be  in  finding  a  proper  analogue  of  the  power  functions  (a)   and  b, 
which  as  ladder  operators  may  not  possess  fractional  powers,  even  though 
they  did  so  when  considered  as  functions  of  a  complex  variable. 

The  Isotropic  Oscillator  in  n  Dimensions 

Using  units  such  that  both  the  mass  and  force  constant  are  unity, 
the  Hamiltonian  for  this  system  is 


H  =  V  (P?  +  x?)/2.  (2.45) 


i  =  l 
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The  eigenfunct ions  of  H  under  the  Poisson  bracket  operation  are 

aj  =  (Pj  -  iXjO/Tz  (2.46) 


with  eigenvalues 


as  wel I  as 


with  eigenvalues 


>^=  i  (2.47) 


a*=  (Pj  +  ix.)//2  (2.48) 


(2.49) 


2 
Hence  one  can  immediately  construct  n   linearly  independent  constants 

of  the  motion  of  the  form  a;aV.  i  =  1,  ..n  and  j  =  1,  ...n,  although 

as  before,  it  is  somewhat  more  convenient  to  deal  with  their  real  and 

imaginary  parts  separately;  in  particular, 

n 
H  =y  aja*  (2.50) 


Dj  =  a,a'f  -  ajaf       i  =  2,  ...n         (2.51) 

K.J  =  a.av  +  a*aj.      j  >  i  (2.52) 

Ljj  =  i[a-faj  -  a.ap        j  >  i  .  (2.53) 

The  latter  three  sets  of  constants  are  closed  under  Poisson  bracket 
and  satisfy  the  following  commutation  relations: 
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im 


->  *VC      .   /-*-       -    *> 


-  (a-fak  -  ajat)6.^  -  (aja,  -  aja*)6.,]     (2.55) 


{Oj.  Dj)  =  0 


(2.56) 


Mata,-Ha.a*)6.,-(ata,.a.ap6,^      (2.57) 

(L.j,  D^)  =  (a^a^  +  a.a*)6,^  -  (aja^  +  aja-^)8j,^ 

-(a|a,  .a.at)5.,  (2-58) 

tDj.  ^kn,^   =  't(^?k  -  ^-"k)^Tn  *  <^K  -  ^j<)^-k 

-  (ata,  -  a,a;)6,,]  .  (2-59) 

Depending  upon  the  values  of  5,  j,  k,  and  m  all  of  the  relations  will 
give  linear  combinations  of  the  functions  Dj,  Kjj .  or  L;j  as  a  result. 
The  relations  also  identify  the  symmetry  group  of  the  problem  [3?]. 
namely  SUp  where  the  n^  -  1  constants  are  considered  as  the  generators 
of  the  group. 


The  Hopf  Happing 

In  general  only  2n-l  independent  constants  of  the  motion  including 
the  Hamiltonian  are  expected  to  be  found  and  not  n2.   Hence  there  should 
exist  n^-(2n-l)  =  (n-1)^  functional  relationships  among  the  n  constants. 
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As  a  consequence,  the  following  set  of  2n-l  constants  have  been  chosen 

for  the  Hopf  mapping; 

the  (n-l)  constants  K,j  where  j  =  2,  ...n, 
the  (n-l)  constants  Lj •  where  j  =  2,  ...n, 

and  the  following  linear  combination  of  H  and  the  Dj 

n 
D=  B,H  +y  B;Dj  (2.60) 

i=2 
where 


6,=^  (2.61) 

I    n 


and 

2 

n 

In  terms  of  the  a;  and  a'j 


B.  =  i  (i  =  2,  ...n)  .  (2.62) 


-I 

i=2 

The  (n-1)^  relationships  between  the  constants  are 


D  =  a^a',  ^  -j-, 
i=2 


(2.63) 


'J        H  +  0 


'-•        H  +  D 


2     2     2 

(H+D)^  -  K,i  -  L,. 
Dj   =  TTTT-^T^ (2.66) 


2(H  +  D) 
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where 

H^'  =  D^  +)   K,f  +)   L,f  (2.67) 


i«2      i=2 


as  may  be  verified  by  substitution  from  the  definitions  given  in  Eqs. 
(2.50)  through  (2.53). 

For  the  Hopf  mapping  let  the  2n  variables  aj  and  a*  be  redefined  as 

a,  =jY^cos  T)e-'P  (2.68) 


a.  =  ^  €.(sin  T)e-'<'j    (j  =  2,  ...n)         (2.69) 

J   -72  J 

where  the  e:  are  direction  cosines  such  that 


n 
j=2 


(2.70) 


Forming   the   ratios 


=   €j(cot   T)e'^^"''J^  (2.71) 


gives  a  complex  projective  space  of  n-1  complex  dimensions.   Regarding 
the  space  as  having  2n-2  real  dimensions  and  performing  an  inverse 
stereographic  projection  a  real  space  of  dimension  2n-l  is  obtained. 
As  in  the  case  of  the  two  dimensional  oscillator  the  sphere  of 
2n-l  dimensions  is  considered  tangent  to  the  2n-2  dimensional  space 
at  its  south  pole.   By  defining  d   as  the  angle  between  a  radius  vector 
from  the  center  of  the  sphere  and  the  diameter  between  the  north  and 
south  poles  2t  can  be  identified  with  0.   The  angles  in  the  planes 


perpendicular  to  the  north-south  diameter  projected  from  the  2n-2  space 

parallel  to  this  line  are  identical  with  the  angles  in  the  space  and 

are  denoted 

«!>!-=  aj  -  p  .  (2.72) 

The  substitution  of  Eqs.  (2,68)  and  (2.69)  in  Eq.  (2,63)  and  into 
the  equations  defining  H,  K) : ,  and  L| :  gives 

H  =  (l/2)\2  (2.73) 

D  =  (l/2)\^  cos  e  =  Z  (2.7^) 

K,j  =  (l/2)\2  sin  eej  cos  Oj  =  Xj           (2.75) 

L,j  =  (l/2)\2  sin  0ej  sin  <Dj  =  Yj          (2.76) 

where  the  K, . ,  L,.,  and  D  have  been  identified  with  the  2n-l  coordi- 
nates in  this  space.   As  in  the  case  for  n  =  2  these  2n-I  constants 
determine  a  point  on  the  surface  of  a  sphere  of  radius  H  in  this  space. 
At  the  same  time  the  remaining  constants  still  define  infinitesimal 
transformations  in  the  respective  coordinate  planes  such  as  rotations 
or  changes  in  the  correlation  or  relative  phases.   It  should  be  noted 
that  only  in  three  dimensions  is  it  possible  to  identify  planes  and 
their  perpendicular  coordinates  in  a  one  to  one  fashion.   In  general 
for  n  dimensions  there  are   (I/2)n(n-l)   coordinate  planes. 


42 


Other  Canonical  Coordinates 

In  finding  canonical  sets  of  coordinates  and  momenta  a  wide  choice 
of  momenta, even  when  the  n  momenta  are  selected  from  amongst  the  con- 
stants of  the  motion, is  available.   The  final  decision  will  depend  upon 
the  particular  application  desired.   However,  there  are  some  simple 
linear  combinations  which  are  often  sufficiently  useful,  to  make  their 
display  worthwhile. 

As  a  first  selection  for  the  canonical  momenta  consider  the  follow- 

n 
ing: 


p,  =  H  =)   ajaf  (2.77) 


i=2 

Pj  =y  ajat  -V  aiaf      (j  -  2.  ...n).        (2.78) 

i=l       i=j 

where  P.  is  the  energy  difference  between  the  first  (j-I)  coordinates 
and  the  last  n-(j-l)  coordinates.  The  coordinates  conjugate  to  these 
momenta  must  then  have  the  form 


(2.79) 


(j  =2,  ...n)  .       (2.80) 

In  terms  of  the  parameters  used  in  the  Hopf  mapping  these  become 

H  =  \2/2  (2.81) 


Ql 

1 

In 

^l^n 
^l^n 

^j 

1 

In 

a*   ,a1 
a.    ,a. 
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P2  =  (\2/2)  cos  d  (2.82) 

j-1  n 

/^2,^^  r   2     \7    2  .  2     V    2  .  2   1 
p.  =  (\  /2)  [cos   T+^   €|Sin  T-\   eisin   tJ 

1=2  l-j 

>  1 
=  (\^/2)  [cos  9  +y   e-sin^  0/2]   (j  =  3.  ...n)  (2.83) 

i»2 

a,  =  l/2(p  -  o„)  (2.84) 

Q2  =  l/2(p  -  02)  (2.85) 

d.   =  l/2(a._,  -  a.)        (j  =  3,  ...n).       (2.86) 

One  might  also  consider  taking  the  total  energy  and  the  energy  in  the 

last  n-(j-l)  coordinates  for  the  canonical  momenta, 

n 

P,  =  H  =y  a. at                    (2.87) 
1=1 


and 


n 

P 

J 


.   =^  ajaf   (j  =  2.  ...n)  (2.88) 


with  the  corresponding  coordinates 


Q:  =  —  In  ^i^^   (j  =  1.  ...n-1)        (2.89) 
J   2,    a^aj^, 


Qn  =  -1"-  (2.90) 

2 1     an 


kk 


which  in  terms  of  the  Hopf  parameters  become 

H  =  >^/2  (2-91) 

n 
p.   =  {\^/2)Y    e?  sin^  T    (j  =  2.  ...n)     (2.92) 

Q,  =  p  -  02  (2.93) 

Qj  =  o.  -  Oj.^i   (j  =  2,  ...n-1)  (2.9^) 

Finally,  one  might  consider  using  the  energy  anci(n-l)  successive 
energy  differences. 


n 

Pi  =  H^  aja; 
1=1 


n 


(2.96) 


P.  =  a.., a*,  -a, at    j  =  (2.  ...n)  .       (2-97) 


The  corresponding  coordinates  are 


1    a)a2  ...a^ 

Ql  =  7-7  ^"  

2ni    a|a2  ■ ■ -^n 

k-1  ^*   n-k+1    n   a-  "^"^ 


(2.98) 


Q,,  =  In  (  ir  (-1)     )  (  t  (4)   )       o  99) 

2n'     j  =  i  ^j        j  =  k  ^j 

k  =  (2.  ...n). 
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In  terms  of  the  Hopf  coordinates 

p,  =  H  =  x2/2  (2.100) 

P2  =  (\2/2)(cos2  T  -  €2  sin^T)  (2.101) 

p.   =   (\2/2)  sin2T(ej.,   -  €j)   (j=3.  ...n)     (2.102) 

n 

Ql  =  (I/n)(p  *Y    o.)  (2.103) 

i»2 

Qj  =  (l/n)[(n-j  +  l)(p^-^  a,^)-(J-0^  ^k^  •       (2.10A) 
k=2         kpj 

In  all  three  cases,  the  transformation  has  been  such  as  to  give 
the  action-angle  variables.   The  effect  has  been  to  map  the  orbit  into 
a  point  on  the  surface  of  a  sphere  of  2n-l  dimensions  whose  radius  is 
the  energy  H.   Since  the  mapping  is  many  to  one,  one  sees  that  as  in 
the  two  dimensional  case  an  entire  orbit  is  mapped  into  a  single  point. 
The  orbital  phase  may  still  be  recovered  by  attaching  a  flag  to  the 
point  representing  the  orbit  which  will  then  rotate  with  constant 
angular  velocity  Qi  ,  since  Q.i  is  conjugate  to  H  and  hence  increases 
1 inear ly  with  t  ime. 

The  second  case  also  has  another  interesting  property,  namely  D2 
is  the  Hamiltonian  for  an  isotropic  oscillator  in  (n-1)  dimensions, 
and  the  other  0;  are  constants  of  its  motion.   Hence  one  may  perform 
another  Hopf  mapping  and  continue  in  this  manner  until  one  has  the  2- 
dimensional  oscillator  again. 
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In  general,  any  of  the  2n-l  independent  constants  from  the  n^-1 

constants  could  be  chosen  as  generators,  and  used  as  coordinates.   It 

should  be  noted  that  the  expressions  for  the  coordinates  could  just 

as  well  be  expressed  in  terms  of  an  arctangent  by  using  the  well  known 

relat  ion 

tan"'  X  =  i  i  In  -L±i  (2.105) 

2      1  +  ix 

which  would  exhibit  the  similarity  between  these  coordinates  and  those 
obtained  by  Goshen  and  Lipkin  [23]. 


The  Anisotropic  Oscillator  in  n  Dimensions 

In  units  where  the  mass  is  unity,  the  Hamiltonian  for  the  n- 
dimensional  anisotropic  harmonic  oscillator  is 


H  =  i£  (P?  +  «^?  x2).  (2.106) 


i  =  l 
where  uj  is  the  frequency  in  the  i   coordinate.   The  eigenvectors  are 
all  of  the  form 

aj  =  (p.    -    i  w.  x.)Af2  (2.107) 

a^  =    (p.   +    i    u,.    Xj)//2  (2.108) 

with  eigenvalues 

\j   =   ±    i    ojj    .  (2.109) 

There  also  exists  n^  constants  similar  to  those  obtained  in  the  iso- 
tropic case,  namely 


^7 


II 


5=1 


*  (2.110) 


Di  =  a..ai  -aj^,  af^,     (i  =  1.  ...n-l)      (2.111) 


(a;  J  a;  '  •»•  a,  J  aj  ')  (2.112) 

Uj-1  Wj-l 

*^2    /o  =.*^  2 


(ajaf)^    (a^ap 

i[a*'^J  a."'  -a"J  af '  ] 

•-I  J  Uj-l  U)j-1 

(a.af)  2   (a^aj)"^ 

Once  again  this  set  is  closed  under  Poisson  bracket  and  satisfies 
the  commutation  relations  of  the  generators  of  SU^. 

It  should  be  noted  that  every  normal  quadratic  Hamiltonian  is 
equivalent  to  some  anisotropic  harmonic  oscillator,  because  they  can 
be  put  in  the  canonical  form  given  in  Chapter  1. 

For  non-normal  Hamiltonians  such  as  the  free  particle 

H  =  -l.P^  .  (2.114) 

2m 

there  is  no  equivalence  with  the  harmonic  oscillator.   However,  symmetry 
groups  still  exist  for  these  cases;  for  example  the  above  expression  is 
invariant  under  rotations.   The  symmetry  of  such  Hamiltonians  has  been 
discussed  in  considerable  detail  by  Moshinsky  [38]  and  F.T.  Smith  [39, 

ko]. 


48 


For  the  Hamlltonlan  given  In  Eq.  (2.114)  (where  m  ■  1)  the  follow- 
ing set  of  nine  quadratic  constants  of  the  motion  exists, 

H  =  (P^  +  P^  +  P^)/2  (2.115) 

X    y    z  \          ■'I 


Ai  -  Pj  -  Pj  (2.116) 

A2  =  Pj  -  Pz  (2.117) 

A3  =  P^Py  (2.118) 

A4  =  P^P^  (2.119) 

A5  =•  PyP^  (2.120) 

Li  =  yPz  -  HPy  (2.121) 

L2  =  zPx  -  xP^  (2.122) 

L3  =  xpy  -  yP^  .                     (2.123) 

The  relations  between  these  constants  are  summarized  In  the  following 
Poisson  bracket  table: 


^9 


H 

Al 

A2 

^3 

A^ 

^ 

h 

•-2 

h 

H 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A| 

0 

0 

0 

0 

0 

0 

2A5 

2A4 

-4A3 

A2 

0 

0 

0 

0 

0 

0 

-4A5 

2A^ 

-2A3 

^^3 
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0 

0 

0 

0 

0 

-A4 

A5 

A] 

^ 

0 

0 

0 

0 

0 

0 

^^3 

-A,-A2 

-A5 

^5 

0 

0 

0 

0 

0 

0 

A2 

-A3 

A4 

L) 

0 

-2A5 

4A5 

Az, 

-^3 

-A2 

0 

•-3 

-L2 

k 

0 

-2A4 

-2A4 

-A5 

A,+A2 

A3 

-4 

0 

Ll 

4 

° 

4A3 

2A3 

-Al 

^ 

-^ 

•-2 

-L) 

0   . 

(2.124) 

These  define  a  group  whfch  is  the  semi-direct  product  [41 ]  of  the  three 
dimensional  rotation  group  and  the  five  dimensional  abel Ian  group  formed 
by  the  five  symmetric  second  rank  tensors.   In  general,  a  free  particle 
in  n  dimensions  will  have  n^-l  linearly  independent  constants  of  which 
n(n-l)/2  will  be  the  generators  of  Rf,.   The  others,  n(n+l)/2  -1  in 
number,  commute  among  themselves  and  are  the  generators  of  the  abel Ian 
group.  These  constants  are  second  rank  tensors  formed  by  taking  pro- 
ducts of  the  momenta. 

The  harmonic  oscillator  Hamlltonlan  will  be  generalized  In  the 
next  chapter  and  consideration  given  to  the  problem  of  a  plane  harmonic 
oscillator  in  a  uniform  magnetic  field. 


CHAPTER  3 

THE  HARMONIC  OSCILLATOR  IN  A  UNIFORM  MAGNETIC  FIELD 

Introduct  ion 

In  the  previous  chapter  it  has  been  shown  how  the  accidental  de- 
generacy of  the  harmonic  oscillator  may  be  accounted  for  by  the  methods 
presented  in  Chapter  1.   Although  historically,  attention  has  most  often 
been  focused  upon  the  accidental  degeneracies  of  the  harmonic  oscillator 
[9]  and  the  Kepler  problem  [k,    5],    the  problem  of  the  cyclotron  motion 
of  a  charged  particle  in  a  uniform  magnetic  field  possesses  certain 
unique  features,  among  which  is  that  It  has  two  linear  constants  of  the 
motion.   It  also  has  accidental  degeneracies  of  its  own  [kl].       Its 

Ham  1 1  ton! an,  _ 

H  =  (P  -  I  A)2/2m.  (3.1) 

where,  for  one  choice  of  gauge 

A=  (-^B^y,  ^iB^x.  0)  .  (3.2) 

is  quadratic. 

As  shown  in  the  previous  chapter,  the  constants  of  the  motion  of 
the  two  dimensional  harmonic  oscillator  K,  L,  and  D  arose  in  a  natural 
fashion  as  products  of  linear  eigenfunct ions  of  the  operator  T^.   The 
eigenvalues  of  Tm  occurred  In  negative  pairs,  so  that  their  products, 
belonging  to  the  sums  of  the  corresponding  eigenvalues,  were  constants 
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of  the  motion.  Moreover,  there  was  a  convenient  mapping,  the  Hopf  mapping, 
which  yielded  a  set  of  canonical  angular  coordinates  <t>  and  ^,  whose  con- 
Jugate  momenta  were  respectively,  H  and  D. 

In  the  present  instance,  one  finds  by  expanding  Eq.  (3.1)  and  using 
the  gauge  in  Eq.  (3-2)  that 

H  -  J.  (p2  ^  p2)  +  li£  (x2  +  y2)  ^  f?o  (yp^  .  ^p^j   (33) 
2m  8mc  2mc 

Effectively  this  is  the  sum  of  a  harmonic  oscillator  Hamiltonian  and  a 
term  proportional  to  the  angular  momentum. 

There  is  one  main  difference  between  the  problem  of  cyclotron  motion 
and  the  harmonic  oscillator.  The  representation  of  the  harmonic  oscil- 
lator Hamiltonian,  applied  to  the  homogeneous  space  H^  '  has  negative 
pairs  of  non-zero  eigenvalues.  All  constants  of  the  motion  are  generated 
by  products  of  the  corresponding  eigenfunct ions ,  for  which  the  eigenvalue 
sum  is  zero.  Thus  the  lowest  order  constants  of  the  motion  are  quadratic, 
belonging  to  h'  ',  and  all  others  are   linear  combinations  of  these  pro- 
ducts. With  the  cyclotron  Hamiltonian  there  are   zero  eigenvalues,  and 
therefore  there  are  elements  of  h'  '  which  are  constants  of  the  motion. 
Whereas  the  Polsson  bracket  of  two  quadratic  constants  is  quadratic,  the 
Poisson  bracket  of  two  linear  constants  Is  a  scalar.   A  rather  different 
kind  of  symmetry  group  results  In  the  two  cases;  for  the  harmonic  oscil- 
lator one  finds  the  unitary  unlmodular  group,  while  for  cyclotron  motion 
one  obtains  a  group  generated  by  the  harmonic  oscillator  ladder  operators. 
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Classical  Zeeman  Effect  for  the  Harmonic  Oscillator 


It  Is  somewhat  more  instructive  to  consider  a  plane  isotropic  har- 
monic oscillator  in  a  uniform  magnetic  field  rather  than  the  cyclotron 
motion  exclusively.   Its  Hamlltonian, 


H  -  i_  (p  -  £  A)^  +  imuj  ^2 


2m 


2m 


(p2  +  p2)  +  El  (0^  +  co2)(x2  ^   y2)  +  ^yp^  .  ^p  ) 


(3.^) 
Is  also  quadratic.   In  Eq.  (3-^)  w  Is  the  Larmor  frequency  eB/2mc  and 
Wq  is  the  natural  frequency  of  the  oscillator.   This  Hamlltonian  reduces, 
In  the  limit  as  cj  -» 0  (B  -»  0)  to  that  of  the  plane  harmonic  oscillator, 
while  as  co^  -» 0  It  reduces  to  that  of  cyclotron  motion. 

The  matrix  representation  of  this  Hamlltonian  considered  as  an 
operator  under  Polsson  bracket  and  denoted  by  T^  Is 


0 

U) 

m( 

u,2 

+ 

"?) 

0 

-to 

0 

0 

m(u)2  +  u)2j 

m 

0 

0 

u 

0 

m 

-U) 

0 

(3.5) 


where  the  basis  is  composed  of  the  monomials  (x,  y,  P  ,  P  ).   The 

X   y 


eigenvalues  and  eigenvectors  of  T^  are 
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Eigenvector 

Eigenvalue 

u 

\[{J  +  u)2j'/2  +  co] 

u* 

-l[(u2  ^  ^)l/2  ^  ^j 

V 

I[(u2  ^.  0,2)1/2   .   ^] 

V* 

-I[(cj2  +  0,2^1/2    .   ^j    ^ 

By  defining 


X  ±  iy 


P  =  Px  ±  iPy  , 
the  eigenvectors  u  and  v  can  be  written  as  follows 


u  -  [m(u)2  +  u2)jl/2  r+  +   I  p+ 


V  =  [m(co2  +  a)2)]l/2  r"  +   I  p' 


(3.6) 

(3.7) 
(3.8) 

(3.9) 
(3.10) 


where  u  and  v*  are  simply  the  complex  conjugates  of  u  and  v  respec- 
tively.  These  four  eigenfunct ions  satisfy  the  following  relation 


(u*.  u)  =  {v*.  v)  =  i+i(u2  ^  Ji^]/2 


(3.11) 


The  constants  of  the  motion  will  be  products  of  eigenfunct ions , 
the  sum  of  whose  eigenvalues  Is  zero.   Hence  one  establishes  the 
quantities  uu'''^,  vv"''',  u'-i^v,  and  u^v*  as  constants  of  the  motion  where 
R  is  a  number  such  that 

R>^|  =  >.2  .  (3.12) 


5^ 


^"'^  2     2  1/2 


In  order  to  display  the  symmetry  group  in  a  convenient  form  the 

following  linear  combinations  are  tatcen  as  the  constants  of  the  motion: 

H  =  [\,  uu*  +  \2  vv*]/i+  yj-   +  J^  (3.15) 

R-1 
K  =  [u^v*  +  u*'^v]/r'/2(uu*)  2  (3.16) 

R-1 
L  =    i[uRv*  -  u*Rv]/r'/2    (iju*)"2~  (3.17) 

D  =    [uu'^  -   Rvv*]/R    ,  (3.18) 

where  H  has  been  written  in  accordance  with  Eq.  (1.29).  These  four 
quantities  satisfy  the  following  relations: 


Tj^(K)  -  T^(L)  =  T^(D)  -  0               (3.19) 

Tk(L)  =  aD  (3.20) 

TjD)  -  aK  (3.21) 

T[,(K)  =  ai,  (3.22) 

a   «  8^u2  +  0)2  .  (3.23) 


where 


The  first  of  these  equations  Is  simply  a  statement  of  the  fact  that  K, 
L,  and  D  are  constants  of  the  motion,  while  the  latter  three  equations 
show  that  the  symmetry  group  of  the  system  is  SU2. 
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If  in  Eq.  (3-12)  R  is  a  rational  number,  then  the  classical  system 
has  bounded  closed  orbits  and  a  quantum  mechanical  analog  exists  for 
the  operators  K  and  L.   However  if  R  happens  to  be  irrational,  then  the 
orbits  are  space  filling. 

Limiting  Cases 

In  taking  the  limit  as  uj  -»0,  i.e.  as  the  magnetic  field  is  turned 
off,  Eq .  (3.^)  becomes  the  Hamiltonian  for  the  two  dimensional  isotropic 
harmonic  oscillator,  the  eigenvectors  become, 

u  «"/m  Uq  r"*"  +  -L-  P*  0.2k) 

V  m 

V  -  NTm  cOq  r-  +  L-  p-  ,  (3.25) 

V  m 

while  the  four  eigenvalues  degenerate  into  two,  namely  ±  IUq,  whence 
R  =  1 .   Expressed  in  terms  of  the  new  u  and  v,  the  four  constants  are 

H  -  (uu*  +  vv*)A  (3.26) 

K  »  uv*  +  u*v  (3.27) 

L  -  i(uv*  -  u*v)  (3.28) 

D  -  uu*  -  vv*  ,  (3.29) 

which  are  the  constants  previously  obtained  for  the  isotropic  oscil- 
lator. The  commutation  rules  for  K,  L,  and  D  also  still  hold. 
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In  considering  the  limit  as  u)q  -»0  one  finds  that 


and 


u  =  NTm  w  r*  +  -U  P+  (3.30) 

V  m 


V  =  -fm   a)r~  +  -L-  P".  (3.31) 


The  eigenvalues  of  u  and  u*  approach  ±  2iw  in  this  limit  while  those 
belonging  to  V  and  v*  both  approach  zero.   Hence,  in  order  to  satisfy 
Eq.  (3.12)  R  must  also  approach  zero.   The  Hamiltonian  in  this  limit 
approaches  that  for  pure  cyclotron  motion  as  in  Eq.  (3-3).   From  the 
values  of  the  eigenvalues  one  immediately  has  two  linear  constants  of 
the  motion,  v  and  v*,  and  one  quadratic  constant  uu*. 

The  constants  can  be  explicitly  derived  by  considering  the  com- 
mutation rules  of  K,  L,  and  D  in  the  limit  as  uJq  -»0. 

Rewriting  the  commutation  relations  explicitly  gives 

,  ruRv*  +  u*^    u^*  -  u'-'^  1   ^  uu^  -  Rvv^'^  ,.    „. 


'uRy*  -  u*K         uu*   -  Rvv*^  „  a     u^*  +  u^'^v  (3-33) 

R-r  (  R-r 


and 


r'/2(,,*)  2        ,     J  R>/2(,,*)  2 


uu*  -  Rvv*   u'^v*  +  u*'^v  I  _  .   ii^v*  -  11*'^ 


R        r'/2(iju  )  2  J       r1/2(uu  )  2 


Multiplying  the  first  of  these  by  R  and  the  latter  two  by  R^^     and  then 
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taking  the  limit  as  cJq  -+0  (R  -♦0),  results  in  the  equations 

I  ((uu*)'/2  (v'Vv),  (uu*)'/2  (v*-v)}  o  8  a)uu*  (3-35) 

I  {(uu*)>/2  (v*-v).  uu*}  =  0  (3-36) 

{uu*.  (uu*)'/2  (v*  +  v)]  -  0  .  (3.37) 

Since  uu   is  simply  twice  the  cyclotron  Hamiltonian,  it  follows  from 
the  last  two  equations  that  both  the  real  and  imaginary  parts  of  v  are 
constants  of  the  motion.   Dividing  Eq.  (3-35)  by  uu*  gives 

^^J1,  1^^   .  8U,  .  (3.38) 

Expanding  Eq.  (3-1)  results  in 

H  =  (P^  +  p5)/2m  +  mJ-    (x^  +  y2)/2  +  u,  (yp   -  xP  )  . 

^   (3.39) 

Thus  the  cyclotron  Hamiltonian  splits  into  two  parts,  one  being  the 
harmonic  oscillator  Hamiltonian  Hq  and  the  second  proportional  to  the 
Z-component  of  the  angular  momentum  L.   Both  of  these  terms  commute 
with  the  total  Hamiltonian  H  =  Hq  +  L-   Hence  another  quadratic  con- 
stant of  the  motion  is 

D  '  Ho  -  L  .  (3.^0) 

For  convenience  in  notation  define 


S  -"^  (v-v*)  .  muy  -  P^  (3.m) 


r 
Q  -^  (v-v*)  -  mux  +  Py  .  (3.^2) 
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With  these  definitions  the  following  commutation  rules  hold: 

(H.  D)  =  {H,  S}  -  (H,  Q}  -  0  (3.^3) 

{S.  Q}  =  2maj  (3./+4) 

{D.  S)  -  2coQ  (3.Zi5) 

(D.  Q)  «=  -2a)S  .  (3.^6) 

These  commutation  relations  coincide  with  those  obtained  by  Johnson 
and  Lippmann  [kl] .      These  authors  have  discussed  the  two  constants  S 
and  Q  in  considerable  detail  and  show  that  they  are  simply  related  to 
the  location  of  the  center  of  the  circular  orbit  and  to  its  diameter, 
which  can  readily  be  seen  In  the  following  manner.   The  canonical 
momentum  expressed  in  terms  of  the  mechanical  momentum  is 

P  =  m V  +  ®  A  .  (3.^7) 

Substituting  for  the  canonical  momenta  in  Eqs.  (3-^1)  and  (3.42)  gives 

S  =  m(2wy  -  v>^)  (3-48) 

Q  -  m(2wx  +  vy)  .  (3.49) 

Evaluating  S  when  Vj^  ™  0  and  Q  when  Vy  =  0  gives  the  center  of  the 

orbit  as 

(xc.  Yc)  "    (Q/2m(J,  S/2mw).  (3-50) 

These  constants  also  determine  the  diameter  of  the  orbit.   Since  the 
orbit  is  a  circle  only  one  of  the  constants  must  be  considered.   For 
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example,  consider  Q.  When  v  takes  on  its  maximum  positive  value,  x 
takes  on  its  minimum  value  and  when  Vy  takes  on  its  maximum  negative 
value,  X  is  a  maximum  and  hence  the  diameter  d  is 

d-  W--min  =  ^^-  (3.51) 

u 

Because  of  the  continuum  of  points  available  for  the  center  for  a 
given  energy  the  degeneracy  of  this  problem  is  infinite. 

Canonical  Coordinates 

As  In  the  case  of  the  plane  harmonic  oscillator  a  set  of  canonical 
coordinates  can  be  found  such  that  the  Hamiltonian  becomes  a  canonical 
momentum.   In  fact,  two  momenta  for  the  problem  are: 

Xi  uu'  +  \o    vv' 

H  =  -' (3-52) 

.,  2    2,1/2  ^^ 


and 


D  =  HHIjlJRvv:  ^  (3.53) 


while  the  coordinates  conjugate  to  these  momenta  are 


Q,  =  \2['n  —  +  >.i  In  ^]Ai  \,\2         (3-5^) 
u  V 


Q2  -      X2[ln  y^  -  \,    In  ^]/l6i    \2    (u)2  +  u,2)'/2        (3.55) 


respect  i vely, 
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A  mapping  similar  to  the  Hopf  mapping  may  also  be  performed  where 

u  =  [\^/2  cos  T  e"'P]  (3.56) 

and 

V  =  [\]/2  sin  T  e-'*']   .  (3.57) 

Under  this  mapping  the  momenta  become 

H  =  X^X^/k   {J-   +  u)^)'^^  (3.58) 

D  =  \]  cos  0  (3-59) 

where 

e  •=  2t  (3.60) 

and  the  coordinates  are 

Q,  =  ^/2^2  (3.61) 

^"''  2     2  1/2 

Q2  =  'J>/8r(u^  +  %)        '  (3.62) 

where 

Y  =  Rp  +  o  (3.63) 

and 

<l>  -  Rp  -  o  .  (3.64) 

Performing  the  mapping  on  the  other  two  constants  of  the  motion  gives 

K  =  X.]  sin  0  cos  (J>  (3-65) 

and 

L  =  \,  sin  0  sin  <P .  (3.66) 

There  also  exists  another  set  of  canonical  coordinates  which  were 
originally  defined  by  Goshen  and  Lipkin  [2k].  Written  In  terms  of  the 
cartesian  coordinates  the  momenta  are  taken  to  be 
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and 


P  =  [(p2  +  p2)/2m  +  m(u)2  +  J)  {k^   +  y^)/2]/{J-   +   u)^)  ' /2(3 .67) 


P0  =  xPy  -  yP^ 
The  corresponding  coordinates  are 


(3.68) 


and 


0=1  tan"' 
2 


q  =  1  tan"' 
2 


[P  P  /m  +  m(u2  +  cj2)  xy] 
X  y 2_ o    

(Px  +  Py)/2m  +  m(uj2  +  '^o)^^^  +  V^^^^ 


(0)2  -^  0,2)1/2  (^p^  ^  yp^)      ^ 
I  (Px  +  P5)/2m  -  m(uj2  +  ui2j(^2  ^  ^2^/2 


(3.69) 


(3.70) 


In  terms  of  these  variables  the  Hamiltonian  has  the  particularly  simple 

form 

H  -  ((^  +  0,2)1/2  p^  .  ^^^  (3.7,) 

from  which  it  follows  that  both  P^  and  Pq  are  constant  in  time. 


Rotating  Coordinates 


The  problem  of  the  plane  harmonic  oscillator  in  a  uniform  magnetic 
field  has  a  certain  uniqueness  when  viewed  from  a  rotating  coordinate 
system.   However,  the  problem  will  first  be  solved  in   plane  polar 
coordinates.  An  extremely  lucid  description  and  tabulation  of  these 
orbits  has  been  given  by  E.R.  Harrison  [^3 ] • 

Assuming  the  direction  of  the  magnetic  field  to  be  in  the  negative 
Z  direction,  the  Lagrangian  is 
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2 

L  =  -  (r2  +  r202)  .  !I!^  r2  -  mwr^Q         (3.72) 

where  the  following  gauge  has  been  chosen, 

Ar  =»  A^  =  0  (3.73) 

and 

The  momenta  conjugate  to  r  and  0  are 


Pr  =  "If-  (3.75) 

and  .  . 

Pq  -  mr  (6  -  u).  (3.76) 


Since  L  Is  cyclic  in  9,    Pq    Is  a  constant  of  the  motion.   The  Hamiltonlan 
Is 


=  (Pj  +  p|/r2)/2m  +  m(u)2  +  wj)  r^/2  +  uPq    .  (3.77) 

and  the  equations  of  motion  are: 

„  2 

Pe       2    2 

p^  =  _£_  .  ni(cj  +  u)o)r  (3.78) 

mr^ 


(3.79) 


r  =  Pr/m  (3.80) 


e    =  _^+  u  .  (3.81) 


In  general,  the  effect  of  imposing  a  uniform  magnetic  field  on  a  system 
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with  a  central  potential  is  to  add  two  terms  to  the  Hamlltonlan,  namely 
a  harmonic  oscillator  potential  which  is  often  neglected  for  small 
fields  [M»],  and  a  term  proportional  to  Pg,  the  angular  momentum. 
The  solutions  for  the  orbit  equations  are: 


_a sin"'    [  br   -2a  i    +  ±_  sin"'    (  ^^   ^   -^  ] 

2  la  I  r2[b2.4aV]'/2  2W  [b^-^aVj^/^ 


a  br^   -2I  "*     .    -1  ^^''o  '^ 

sin"'    {  -^ — s 701/9    ) sin"'    {  — 5 5-5    \  Jn    ^ 

2|«|  rJ[b2.4aV]'/2  2W  [b^-UaVl'^^ 

(3.82) 

where  2    «  « 

b  =  £_+  w2r2  +  r2  (3-83) 

a  =  P0/m  (3.8^) 

W=  (0,2  +  0,2)1/2  .  (3.85) 

The  subscripts  on  the  coordinates  and  velocities  denote  initial 
values.   These  orbits  are  plotted  in  Fig. 3-1  through  3-^.  I"  al '  cases 
units  have  been  chosen  such  that  m  <»  q  =  c  =  1.   Fig.  3-1  and  3.2 
show  the  high  and  low  field  orbits  respectively,  for  the  same  set  of 
initial  conditions  and  Wq.   Fig.  3-3  and  3-^  show  the  orbits  for  a 
fixed  magnetic  field,  but  for  different  values  of  the  initial  tan- 
gential velocity.   In  general  the  orbit  will  be  a  processing  ellipse 
for  Pq  ^  0   and  will  be  a  hypotrachoid  similar  to  Fig.  3-1  for  Pq  <  0. 
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Parameters 


Initial  Conditions 


u)  =  63A 


10   «  4 

o 


e  =  0 

o 


r  ■  0 
o 


Fig.  3-1   Particle  orbit  in  a  strong  uniform  magnetic  field 
and  harmonic  oscillator  potential. 
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In  a  rotating  coordinate  system  defined  by 

e  =  e  -  cot  (3.86) 

r  -  r  ,  (3.87) 

where  w  is  the  Larmor  frequency,  the  Lagrangian  is 

The  canonical  momenta  in  the  rotating  coordinates  system  are: 

P-   =  mr  (3.89) 

P-  =  m?^5   .  (3.90) 

o 

L  is  also  cyclic  in  0  and  hence  P-  is  a  constant.   The  Hamiltonian  in 

0 

the  rotating  coordinate  system  is 


H  =  P^r*  P-  0  -  L 


2       „  2.-2.  ,,     ,2  2.-2 

+  PS  /r  )/2m  +  m(u)  +  o)  ) 
r     o  ^      o 


(P;  +  Pq  /?  )/2m  +  m(u  +  <)r  /2         (3-91) 


and  hence  the  equations  of  motion  are 


p.  =  _L  -  m(cj2  +  cj2)p  (3.92) 

r  _9  u 

Pg  =  0  (3.93) 


?    =  P?/m  (3.9M 


e ^     .  (3.95) 

mp2 
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which  are  the  equations  of  motion  for  a  plane  Isotropic  harmonic  oscil- 
lator with  force  constant  m(u)2  +  u^)  .  In  general  the  transformation  of 
a  Hamiltonian  with  a  central  force  from  a  stationary  coordinate  system 
to  a  rotating  coordinate  system  will  subtract  a  term  proportional  to  Pq. 
Generalizing  then,  it  can  be  stated  that  a  Hamiltonian  with  a  central 
potential  and  with  a  uniform  magnetic  field  present,  when  viewed  from 
a  rotating  coordinate  system,  has  a  harmonic  oscillator  potential  added, 
that  is 

H  =  f  +  V(r)  +  1  mu)^  r^  (3-96) 

where  T  Is  the  kinetic  energy  in  the  rotating  coordinate  system  and  u 
Is  the  rotation  frequency  of  the  coordinate  system  with  respect  to  the 
fixed  system.   In  the  present  case  since  \/(r)  is  a  harmonic  oscillator 
potential  one  simply  has  a  harmonic  oscillator  with  a  larger  force  con- 
stant. 

The  solution  for  the  orbit  defined  by  the  Eqs.  (3.92)  through 
(3.95)  is 

-   .     a  br^-2a^         a       ,      br^-2a^ 

e-eo'-—   sin-'  {   ,  ,   ,  ^  ^    } sin-'  [   ,   °  !  .  ...  ]. 

2|a|         r^Cb^-^aV]'/^    2|a|         r^[b^-^aV] '^^ 

where  W  Is  defined  in  Eq.  (3.85)  and 

a  =  Pg/m  (3.98) 

b  =  a2/r2  +  W^  r^  +  f2  ^3  g^j 

Equation  (3-97)  Is  the  equation  for  an  ellipse,  regardless  of  the 
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magnitude  of  o).  The  orbits  corresponding  to  Fig.  3.1  and  3.2  in  a 
rotating  coordinate  system  are  shown  in  Fig.  3-5  and  3-6. 

Discussion  of  the  Constants  of  the  Motion 

Of  the  four  constants  of  the  motion  defined  in  Eqs.  (3-15)  through 
(3.18)  one  is  H,  the  energy  of  the  system.   Another  constant  is  the 
angular  momentum  Pq   which  is  a  linear  combination  of  H  and  D 

Pe^  JiL-H  -     '  -^  ^     D.  (3.100) 

Since  only  three  of  the  constants  are  independent  one  would  like  to 
discover  one  other  constant,  independent  of  the  two  above,  which  has 
a  physical  meaning.   This  can  be  done  by  considering  the  constants  of 
the  motion  in  the  rotating  coordinate  system  and  then  transforming  back 
to  the  original  coordinate  system. 

In  the  rotating  coordinate  system  two  constants  of  the  motion  are 
aa*  and  bb'"',  where 

a  =  Pjj  -  \{J-   +  u)2)'/2  -  (3.101) 

and 


b  -  P-  -  1(0)2  ^  ^)l/2  -^  (3.102) 


and  where  the  mass  is  taken  to  be  unity.   By  a  simple  rotation  the  axis 
can  be  oriented  so  that  the  semi-major  axis  of  the  ellipse  lies  along 

X  and  the  seml-mlnor  axis  lies  along  y.   In  this  case  then  aa'  is  pro- 

-2 

portional  to  the  maximum  value  which  x  attains,  because  since 
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o 
I 

CD 
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aa*  =  P-2  4-  (uj2  +  cj2)  ^2  (3.103) 

is  a  constant,  it  can  be  evaluated  when  P.  =  0.   Similarly  bb*  Is  pro- 
portional  to  y  ^      .      Alternatively,  one  can  say  that  aa  and  bb''  give 
the  boundaries  of  the  orbit  i.e.  aa'  gives  the  maximum  radial  distance 
the  particle  can  attain  and  bb*  is  the  minimum  radial  distance. 

These  two  radii  are  also  constants  of  the  motion  in  the  stationary 
system  since  the  transformation  only  involved  the  angle.   However,  it 
can  explicitly  be  shown  that  these  are  constants. 

By  transforming  from  the  rotating  coordinate  system  to  the  station- 
ary one  aa'''  can  be  written  in  terms  of  the  eigenvectors  and  their  con- 
jugates as 

aa"'^  =(l/i+)(uu'"'  +  vv"'^  +  uv'^  e^''^^  +  u*v  e'^'^^tj  ^   (3.104) 

In  the  stationary  system  the  time  rate  of  change  is 

^  (aa*)  -  (aa*,  h}  +  ^  (aa*) 

-(lA)[-2iaxjvV'^t  +  2iu)u*ve-2''*'t 

+2iu)uvV'^t  -  2iuju*ve"2''^t]  »  0,    (3.105) 

and  hence  aa*  is  a  constant  of  the  motion  in  the  stationary  system  also. 
A  similar  proof  also  shows  bb*  to  be  a  constant.   However  aa*  and  bb* 
are  not  independent,  but  are  related  to  one  another  through  H  and  D. 

Even  though  a  good  physical  interpretation  cannot  be  given  to  the 
constants  K,  L,  and  D  their  effect  on  the  orbit  under  Poisson  bracket 
can  be  calculated  i.e.,  the  infinitesimal  change  each  produces. 


7^ 


Before  calculating  the  effects  of  the  constants  on  the  orbit,  it 
is  desirable  to  write  the  orbit  equation  in  terms  of  the  constants  of 
the  motion  and  the  time.   This  can  be  done  by  defining  a  complex  vector 

7  =  X  +  iy  (3-106) 

which  can  be  written  in  terms  of  the  eigenvectors  by  inverting  Eqs. 

(3.9)  and  (3.10) 

r  »  (u  +  v'0/2c  ,  (3.107) 

where 

c  =  [m(u)2  +  u^)]'/2.  (3.108) 

Employing  the  polar  forms  of  both  the  eigenvectors  in  Eqs.  (3.56)  and 
(3.57)  and  of  the  constants  In  Eqs.  (3-59),  (3.65)  and  (3-66),  u  and  v* 
may  be  expressed  in  terms  of  the  constants  as 

u  =  (R/2)'/2  [(k2  +  l2  +  0^)1/2  ^  D]'/2  e-'P      (3.109) 
v*«.  (R/2)'/2  [(k2  ^   L^  +  d2)1/2  .  q]W2   ^\a  ^3,10) 

However  from  Eqs.  (3-63)  and  (3.64) 

p  =  (1/2R)(Y  +  <D)  (3.111) 

and 

o  =  (l/2)(f  -  <t>)    .  (3.112) 

Hence 

u  -  (R/2)'/2  [(k2  ^   l2  ^  d2)'/2  +  D]'/^  e-'*/2R  e-'^/^R 

(3.113) 
and 

v*.(Fl/2)'^'  Uk'  *  l'  *  0')'''   .  01'''  e-'"^  e""^    . 

(3.114) 
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From  Eqs.  (3-65)  and  (3-66)  It  Is  apparent  that 

<D  =  tan"'  i  .  (3.115) 

K 

The  time  dependence  Is  in  ^  since  Y/2\2  'n  Eq.  (3.61)  is  the  coordinate 
canonical  to  H  and  hence  from 

d/dt  (^/2\2)  =  {W2\2.  H)  -  I  (3.116) 

it  fol lows  that 

^=  2\2(t  -  t^)  (3.117) 

or 

^  -  2R\,(t  -  t^)  .  (3.118) 

Substituting  Eqs.  (3.115).  (3.117),  and  (3.118)  into  Eqs.  (3.113)  and 
(3.11^)  and  then  substituting  these  into  Eq .  (3.107)gives 

7=  i_  (R)1/2{[(k2  +  l2  >  D2)l/2+D]'/2e(-l/2R)tan-l  (L/K)  e'i^Kt-to) 
2c  2 

+  [(K^  +  L^  +  D2)'/2.p^]/2^(-i/2)tan-l  (L/K)  ^i^2(t-to), 

=  '^-  +  ^  .  (3.119) 

where 

2c  2 

(3.120) 


7_  .  1_  (a)'/2[(K2  +  L^  ^  D2)'/^D]'/2e(-'/2R)tan-'  (L/K)  ^-U,(t-to) 


and 

7^  =  i_  (R)'/2[(k2  +  L^  +  D')'/2.o]'/23(-'/2)tan-l  (L/K)   J^2(t-to). 
2c  2 

(3.121) 
Equation  (3.119)  expresses  the  orbit  equation  as  the  sum  of  two  contra 
rotating  vectors,  whose  frequencies  of  rotation  are  \j  and  Xj,  and 
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whose  amplitudes  are  functions  of  the  constants  of  the  motion. 

The  Poisson  bracket  of  the  constants  K,  L,  and  D  with  7  are: 


(K.  7)  =  ^ 
2 


IKD 


[(K^  +  L^  +  d2)'/2.d]   k^  +  L^ 


iKD 


[(K^  +  L^  +  D^)'''^D]    R(K^  +  L^) 


(3.122) 


(L,  r) 


a 

+  — 


r/  2    2    2J/2   ,     ,2    2, 
[(K  +  L  +  D  )    +D]    R(K  -^  L  )_ 


(3.123) 


(D.  r) 


\a  _  Aa  _ 
2R  '■-'V  '■■» 


(3.12^) 


where  a  is  a  constant  defined  in  Eq.  (3-23). 

The  infinitesimal  change  induced  by  the  constants  is  that  each 
changes  both  the  amplitude  and  phase  of  r  and  r_,  while  preserving 
the  sum  of  the  squares  of  their  sum  and  difference.   If  the  orbit  were 
an  ellipse,  this  last  statement  would  be  equivalent  to  stating  that  the 
sum  of  the  squares  of  the  semi  axes  of  the  ellipse  is  a  constant.   The 
proof  of  the  statement  goes  as  follows: 

Before  the  infinitesimal  change  the  quantity  is 


Ir^+  r.|2  +  |r^  -  rj^  =  2(|rj2  +  |rj2)      (3.125) 
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and  after  the  change    it    is 

|r^(l   +  €)   +   r.(l    +  5)|2  +   |r^(]    +  e)    -   r.(l    -  5)|2 
-   2[|r    I^  +    |rj^   +   (€  +  €*)    lr_^|2  +   (6  +  d'")    |rj^]      (3.126) 

where  second  order  terms  have  been  ignored.   The  difference  between  the 

two  terms  is  ^      o         u      o 

2[(€  +  e')    |r^|2  +  (5  +  5*)  irj2]  . 

which  is  zero  in  all  three  cases.   For  example,  if  the  infinitesimal 
change  is  induced  by  K,  one  has  for  the  difference 


2[(e  +  6--^)|r^|2  +  (6  +  6*)|rj2] 


2a 
2 


2L  frj' 


2L  |r.r 


[(k2.l2.d2)'/2.d]   [(k2.l2.d2)'/^d: 


„  2aR 
"  8c 

-  0. 


[(K^  +  L^  +  d2)'/2.d]     [(^2  +  L^  +  d2)'/2+oJ 


(3.127) 


Gauge  Transformations  in  Uniform  Magnetic  Fields 


In  uniform  magnetic  fields  and  in  the  absence  of  other  external 
potentials  the  following  fact  concerning  a  charged  particle  in  this 
field  is  true:   The  change  in  momentum  in  going  from  one  point  to 
another  Is  Independent  of  the  path  taken  between  the  points.   The 
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proof  is  as  follows:   Let  P  denote  the  canonical  momentum,  ir  the  mechani 
cal  momentum,  B  the  uniform  field,  and  A  the  vector  potential  such  that 


and  since  B  is  uniform 


Since 


B  "  Vx  A  (3.128) 


A  -  i  B  X  r   .  (3.129) 


^^  B.-7  xB    .  (3.130) 

at   c 


and 


one  has 


dt 

/   dn  =  J  /   (f  X  B)  dt 

^2  _   _ 

=  ^  /  d?  X  B 
c  •' 

^2 


c 
r 


Hence 


and 


c 


(3.131) 


^/  dA  .  (3.132) 

1 


Art  "  -  ^  AA  (3.133) 

c 

AP  -  -  ^  AA  .  (3.13M 


Equation  (3.134)  states  that  the  difference  in  momenta  between  any  two 
points  is  proportional  to  the  difference  in  the  vector  potential  evalu- 
ated at  these  points.   From  this  it  is  inferred  that  the  choice  of  a 
certain  gauge  is  equivalent  to  picking  the  zero  of  momentum. 
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One  can  also  see  why  a  translation  of  coordinates  must  be  accompanied 
by  a  gauge  transformation.  When  the  translation  occurs  the  zero  of  momen- 
tum changes  and  hence  this  change  must  be  subtracted  from  A  to  give  the 
the  same  zero  of  momentum.  However,  the  change  in  A  is  simply  a  gauge 
transformat  ion. 

In  order  to  show  that  a  gauge  transformation  can  induce  a  translation, 

first  consider  the  case  of  a  free  particle  in  a  uniform  field  with  the 


following  gauge         A-^(-y  x  0) 

2 


(3.135) 


In  units  where  m  and  co  -  eBQ/2mc  are  one,  the  constants  of  the  motion  are 

Px  -  y  =  «  (3.136) 

and 

Py  +  ^  =  8  •  (3.137) 

Using  the  fact  that 

P-^+fA.  (3.138) 

the  location  of  the  center  of  the  orbit  is  found  to  be 

^\'   Vc)  "   (f.  -f)  •  (3.139) 

However  with  the  gauge 

D 

A  -  -2  (  -(\  +  y).  ^i  +  X,  0).  (3.1^0) 

even  though  the  two  constants  are  the  same,  one  finds  the  center  of  the 
orbit  to  be 

("c  Vc)-(2^  ^)  .         (3.141) 
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Hence  the  center  of  the  orbit  has  been  translated  -g,/2  units  In  the  x 
direction  and  \/2  units  in  y. 

The  center  of  the  cyclotron  orbit  may  also  be  rotated,  however 
this  Is  not  accomplished  by  a  simple  gauge  transformation.   Instead 
one  must  rotate  the  vector  A  and  also  rotate  the  coordinate  system. 
For  the  vector  potential  defined  in  Eq.  (3-135)  the  center  of  the  orbit 
is  given  In  Eq.  (3.139).   A  rotation  of  the  vector  field  defined  by  A 
g  Ives 


h-'f 


cos  e  sin  e  0 

-sin  0   cos  d        0 
0        0     1 


-2.  {    -y  cos  9  +  X  sin  d,      y  sin  0  +  x  cos  6,  0).   (3.1'*2) 


The  rotation  of  the  coordinates  defined  by 


and 


X  ■  X '  cos  0  -  y '  s  in  0 


y"x'  sin0+y'  cos  0 


(3.143) 
(3.1M*) 


transforms  Ad  into  A'  where 


A-  »  _°  (-y',  X',  0)  . 


(3.1^5) 


Hence,  the  center  of  the  orbit  Is 


and 


(x' ,  y ')  ■  ( — , ) 

^  c   c     2     2 


S  „  |1  cos  0  +  |1  sin  0 
2   2         2 


(3.146) 
(3.147) 
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1=  .  |1  sin  e  +  |1  cos  0,  (3.1^8) 


which  shows  that  the  center  of  the  orbit  has  been  rotated  with  respect 
to  the  original  coordinate  system,  while  the  diameter  of  the  orbit  has 
remained  unchanged. 

The  vector  field  can  also  be  rotated  without  rotating  the  coordi- 
nates and  still  preserve  the  curl  If  at  the  same  time  the  vector  is 
dilated  by  a  factor  of  sec  0.   The  new  vector  potential  is 

D 

A  =  —  sec  9   (-y  cos  0+xsin0,   ysin0+x  cos  d,   0). 

2  (3.149) 

The  transformation  leaves  the  center  of  the  orbit  invariant  and  in  fact, 
Is  equivalent  to  a  gauge  transformation  where  the  vector  A^ 


A.  =  ^  tan  e(x.  y.  0).  (3.150) 

whose  curl  is  zero,  has  been  added  to  the  original  vector  potential  A. 


CHAPTER  k 

THE  KEPLER  PROBLEM 

The  Two  Dimensional  Case  in  Polar  Coordinates 

Tiie  Ham  i  I  ton  i  an  for  tlie  two  dimensional  Kepler  problem  in  polar 

coordinates  is 

H  «  (1/2)(P;  *   p;/r2)  .  1/r  {k.\) 

where  the  mass  and  force  constant  have  been  chosen  equal  to  unity. 
The  corresponding  Hami 1 ton-Jacobi  equation  is 

(^,^  .  '    (^)^  -  2(W  .  1)  =  0.  (^-2) 

br  r^        ^  r 

where  W  is  the  total  energy  and 

S(r,  <f)  =  Sr(r)  +  S^(<t>)    .  (4.3) 


Hence 

p  «^  (^•^) 


and 


'*  '  i  ("5) 


Since  the  partial  differential  equation  (^+.2)  has  no  <J>  dependence 
bS/^   "  P(K  Is  a  constant  and  is,  in  fact,  the  angular  momentum  of  the 
system.   Solving  Eq.  (4-2)  for  3s/3r  »  P^.  one  finds 

|S  „  p    [2W+  2  -  li]  (4.6) 

or  r   r'^      ■ 
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The  action  variables    J^  and   J^  are  given   by   tlie    integrals 

Jr  =  55  P^dr  (4.7) 

J^'  P  P4,d<t»    .  (4.8) 

Equation  (4.8)  can  be  integrated  immediately  since  the  integral  over 
one  cycle  is  simply  an  integral  from  0  to  2n,    hence 

J^  =  2«P^  .  (4.9) 

The  integral  in  Eq .  (4.7)  can  be  evaluated  by  a  contour  integration 
as  described  in  Born  [35]  and  the  result  Is 

2n 

The  angle  variables  Wp  and  w^  which  are  canonical  coordinates  and  con- 
jugate to  J^  and  J(j,  are  given  by 


and 


Upon  performing  the  Integrations 

w^  =   (u   -  e  sin  u)/2«  (4.13) 

and 

W(j  -    [u   -   6   sin  u   -   tan"'(y/x)    +  <t>]/2jr,  (4.14) 


m 


where  the  elliptic  anomaly  u  Is  defined  through 

r  =  a(l  -  €  cos  u).  ('♦•15) 

a  being  the  semi-major  axis  and  €  being  the  eccentricity  of  the  ellipse, 


a  «  -  W/2  (4.16) 


€^  =.  1  +  2wpJ   .  (4.17) 


By  choosing  <t>      to  be  the  azimuth  of  the  perihelion, 


X  =  r  cos  (4"  -  ^^)  (4.18) 

y  =  r  sin  (<l>  -  *)•  )  ,  (4.  19) 

and  Eq.  (4.14)  reduces  to 

w.  =  (u  -  €  sin  u  +  <(>  )/2jt  .  (4.20) 

(p  o 

Solving  Eq.  (4.10)  for  the  energy  in  terms  of  J  and  J^ 

,  2 
W=  ^ (4.21) 

(Jr  +  V 
and  hence  the  frequencies  v  and  v^  are   degenerate  i.e., 


r     <*> 


^=  ^  =  ^=  ^  •  ^"^""^ 


Thus  the  quantities  aa' ,  bb  ,  ab' ,  and  a  b  are  constants  of  the  motion 

where  J   1/2  2jtiw 

a  =  (2^    e    '  (4.23) 
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and  b^A'/'e^"'^.  (4.2^ 

2  It 

These  quantities  satisfy  the  following  Poisson  bracket  relations 

{a,  a*)  =  (b,  b*}  -  I  (^.25) 

where  one  uses  the  fact  that  the  J's  and  w's  form  a  set  of  canonical 

coordinates  and  momenta. 

Taking  the  real  and  imaginary  parts  of  these  constants  and  the  sum 

and  difference  of  the  real  ones  gives  the  following  four  constants  of 

the  mot  ion, 

H  »  aa*  +  bb"'^  (^.26) 

D  =  aa'''  -  bb''^  (^.27) 

K  =  ab*  +  a'*b  (4.28) 

L  -  f(ab*  +  a*b).  (4.29) 

Of  these  constants  H  commutes  with  the  other  three  and  is  a  function  of 
the  Hamiltonian,  namely,  it  is  inversely  proportional  to  the  square 
root  of  the  energy.   The  three  constants  K,  D,  and  L  satisfy  the  com- 
tation  relations  of  the  generators  of  the  group  SUo  or  R,. 


mu 


(K,  0}  -  2L  (4.30) 

(D,  L)  "  2K  (4.31) 

(L.  K]  =  20  .  (/4.32) 

Although  they  have  no  simple  interpretation  as  constants  of  the  motion, 
one  can  make  a  mapping  of  these  constants  which  preserves  the  commuta- 
tion rules.   The  new  constants  can  be  interpreted  as  the  angular 
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momentum  and  the  two  components  of  the  Runge  vector  [11].   The  mapping 

(^.33) 


K  -.K'  =  K(2  +  H±D)l/2 
H-D 


L  .^L'  -  L(2  +  H±D)l/2 
H-D 


(k.3k) 


D  -»D'  =  H  -  D  . 


These  constants  satisfy  the  rules 


{K',  D'}  =  -2L' 
(D',  L')  =  -2K' 
(L',  K'}  =  -2D' 


(^.35) 


(4.36) 
(^.37) 
(^.38) 


Otherwise  they  are  a  direct  consequence  of  the  commutation  rules  for 

K,  L,  D,  so  that  this  mapping  is  valid  for  any  SU2  group  where  use 
must  be  made  of  the  identity 

h2  =  K^  ■»■  L^  +  d2  .  (4.39) 

K'  Is  essentially  the  x  component  of  the  vector  P  where 


P  =  RAT^   =  (-cATTi) 


cos  <)>, 


s  in  <t>. 


(i+.40) 


and  R  is  the  Runge  vector,  while  L'  Is  the  y  component  and  D'  is  the 
angular  momentum. 
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The  Kepler  Problem  In  Three  Dimensions 

Using  the  same  units  as  in  the  previous  case,  the  Hamiltonian  for 
the  Kepler  problem  in  three  dimensions  Is 

2      ? 

'   2   Pe    P*      I 
"--  (Pr^-7%  \     )  --  •      (^.^0 
2       r^   r^  sin2  e    r 

Setting  up  the  Hami I ton-Jacobi  equation  and  separating  gives 


as    '  ' 


ae  V  ^-    sln2  6/ 


and 

as        2   ai 


P^  »  -.=  (2W  + ^)'/2^  (1^  i^) 


where  Oq  and  o^  are  separation  constants,  a^   being  the  magnitude  of 
the  Z  component  of  the  angular  momentum  and  Oq   the  magnitude  of  the 
total  angular  momentum.   Since  the  potential  is  central  the  orbit  Is 
planar  and  the  normal  to  the  plane  Is  necessarily  parallel  to  the 
angular  momentum.  Denoting  by  B  the  dihedral  angle  which  the  orbital 
plane  makes  with  the  x-y  plane  one  can  deduce  that 

cos  B  =  ^  .  (4.^5) 

Solving  for  the  action  variables  gives 

J4,  «.  2na^  -  2nP^  {k.kS) 
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Jq  -  2(00  -  a^) 
2 


-2W 


J0  -  J4, 


(k.kQ) 


and  from  Eq.  (4.48)  It  follows  that 


-2^ 


(j^  +  Je  +  J({,)' 


(4.49) 


In  the  plane  of  the  orbit  rectangular  coordinates  M-  and  v  can  be 
introduced,  where  the  \l  axis  is  the  major  axis  of  the  ellipse  and  the 
origin  coincides  with  the  center  of  force  and  the  origin  of  the  polar 
system.   In  this  system 

H  =  r  cos  (Y  -  Yq)  (4.50) 


V  =  r  sin  (^  -  ^o), 


(4.51) 


where  ^  is  an  angle  measured  from  the  M-  axis  to  the  particle  and  ¥q 
Is  the  angle  to  the  perihelion. 
The  angle  variables  are 


3Pr 
w_  =>  /  — -  dr  =  (u  -  e  sin  u)/2n 


(u  -  e  sin  u  -  tan-1  y.  +  sin"'  S2^)/2n 
^  [I  sin  B 


(4.52) 


(4.53) 
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-  (u  -  e  sin  u  -  tan-'  ^^  +  sin"'  S£L^   .  sin"'  ^^i-g  +  0)/2« 
^  M-        sin  B        tan  B    ' 

(^.54) 
where  u  Is  defined  as  In  Eq.  (^.15)  and  €  as  In  Eq.  (4.17).   Using  Eqs. 
(4.50)  and  (4.51)  it  follows  that 

tan-^  V  .  ^j,  .  Y  (4  55J 

It  follows  from  Fig.  4.1   that  sin"'  S°l-2.   is  the  projection  of  the 

tan  B 

angle  H   onto  the  plane  in  which  <t>  is  measured  i.e.,  the  x-y  plane,  which 

shall  be  denoted^  .   It  also  follows  that  sin"'  £2L-£  Is  the  angle  1, 

sin  B 

where  the  line  of  nodes  has  been  chosen  to  coincide  with  the  major  axis  of 
the  ellipse.   Hence,  the  line  of  nodes  makes  an  angle  (4>  - -^)  with  the 
X-axis.   Rewriting  Eqs.  (4.53)  and  (4.54)  results  in  the  following  set 
of  equations, 

v^  =  (u  -  e  sin  u  +  ^o)/2n  (4.56) 

W4,  =  (u  -  e  sin  u  +  ^Q  +  <^  -  ^)/2jt  .        (4.57) 
By  defining 

a  =  (-L)    e    "^  (4.58) 

2it 


b  »  (:i^)'/2  e2«Iv^  (4.59) 


In 


J*    1/2     27tIwij,  ^,    ^   ^ 

c  =   (-i)  e       ^  (4.60) 

In 
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Fig.  ^.1   Three  dimensional  coordinate  system  for 
the  Kepler  problem. 


91 


one  finds  the  following  nine  constants  of  the  motion:   aa*,  bb*,  cc*, 
ab  ,  a  b,  be  ,  b  c,  ac  ,  and  a'c.   These  constants  are  the  generators 
of  the  group  SU^ .   The  constant  angles  B,  ^^  and  <>  - -^  and  also  the  sum 
♦  - -^  +  Hq   can  be  expressed  in  terms  of  the  constants, 


sin  ^o  -  iii.^— ^-i-M.  (4.61) 

°   2(aa*bb*)'/2 


^     ab  4-  a  b  ,,  ^„, 

cos  ^o fTf  (^.62) 

2(aa*bb*) 


slnC-^.  '\'^l\f,%  (^.63) 

2(bb'*cc''')'/2 

cos(»  --^)  =  ^^\\^y,  (^.64) 

2(bb*cc^')'/2 


in(^  --+YJ  =  '(^"^^  -  ^y,l 
°  2(aa*cc*)'/2 


(^♦•65) 


and 


cos(..^.,„).^— ^-,,  C-^^ 


COS  B  -  --£S ^  .  (i,.67) 

bb*  +  cc* 


Generally,  the  rotation  group  in  four  dimensions  R^,  is  tal<en  to 
be  the  symmetry  group  of  the  Kepler  problem.   Its  generators  are  the 
three  components  of  the  angular  momentum  vector  L  and  the  three  compo- 
nents of  the  normalized  Runge  vector  P  which  is  defined  in  terms  of  R 
in  Eq.  {k.kO).      The  rectangular  components  of  P  are 
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•2W 


For  L  one  has 


cos  ^o  cos  (♦ 
-cos  ^Q  sin  (♦ 
sin  B  sin  ^„ 


"-)    -   cos  S  s  i  n  ^Q  sin  (<t>  -  -^) 
"-)    -   cos  S  sin  ^Q  cos  (<t>  -  ^) 


(4.68) 


In 


sin  B  sin  (<t>  -  -^) 
s  in  S  cos  (♦  -  -^) 
cos  B 


(4.69) 


Using   Eqs.    (4.61)    through   (4.67)    one  can  express   the  generators  of  R/^ 
In   terms  of  the  generators   of  SU3   as   follows: 

^'"!^f?Vy^'^!     I(ab*  -.  a*b)(bc*  *  b*c)(bb*  +  cc*) 


""       4bb*(cc*)'/2(bb''  +  cc*) 


+  (ab*  -  a*b)(bc*  -  b*c)    cc*] 


(4.70) 


p     .    -l(aa*-H2bb*-.2cc-'^)^^^-^3t^.v  ^  a*b)(bc*  -   b*c)(bb*  +  cc''') 
y       4bb-*(cc*)'/2(bb*+  cc*) 


+  (ab*  -  a*b)(bc*  +  b*c)    cc*] 


(4.71) 


1  /o 
i(aa*  >  2bb*  +  2cc*)  ,,,-k   .    „     •>>  1 /o,    .*         --v.v  ,,    ,,v 

P     »  -i — — T ;- — ^ (bb     +2cc)'/'^(ab     -  a  b)  (4.72) 

^  2(bb*  +  cc*") 


i(bc*  -   b*c)(bb*  *  2cc*)^^^ 
2(cc*)'/2 


(4.73) 


_   (bc*-^b-'-c)(bb-^-H2cc'^)'/2 
•-y  2(cc*)'/2 


(4.74) 
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and                                                   L^  -  cc*  .  (4.75) 
These  six  quantities   satisfy   the  relations 

^^l'  Pj5  -  -  ^ijk  Pk  ('+•76) 

a,.    Lj}  -  -  €,.^  L,  (ii.77) 

^Pi'   Pj^  -  -  ^ijk  k  (^.78) 

which  define  R^,  once  again  directly  from  the  commutation  rules  for  SU3, 

The  process  can  be  inverted  and  the  generators  of  SU3  can  be  ex- 
pressed in  terms  of  the  components  of  P  and  L.  These  are 

..*.^.(,2,,V/2.,^  (-..79) 


-N/-2W 

bb'^ML^.L^)'/^  (4.80) 

cc*  -  L^  (4.81) 


'*  "    IT^  "   ^'■'  *   '?^"'-   ^^^"'<^x  -   l2)>/2(cos  1^   -    i    sin  ^ 


■2W 

(4.82) 
be*-    tl-^d-^  +   l2)'/2]1/2    (^^5    (^   .^)    .    ,    sin    (4,    ._))  (1^83) 

ac*-    [7=i=-(L2  +   l2)'/2.l^]1/2lJ/2[,,3(,   .^.,^^).i    sin(*   --*  ^^ 


where 


2w 

(i*.84) 


"-x 


^  -  ^'  tan"'  -^  (4.85) 
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^Q  -  tan 


P^LydJ  +  d  ^   l2)'/2 


•1  ^z'-y 

P  (L^  +  L^)  +  P  L  L 
■^x^  X   y'    z  X  z 


and 


S  =  tan 


^2  ^  ,2J/2 
-1  (Lx  +  Ly) 


(^♦.86) 


(4.87) 


The  Two  Dimensional  Kepler  Problem  in  Parabolic  Coordinates 


In  parabolic  coordinates,  defined  by 


X  a  ^i  V 


y  =  (n^  -  v2)/2 


and  where  m  =»  e  =  1,  the  Kepler  Hamiltonian  is 

Pu  +  P^  -  4 


H 


[i     ^  ry 


2(^2  +  v2) 


From  the  Hami 1  ton- Jacob i  equation 


where 


P,  =  (2Vii2  +  2a,)  ^/2 


P  =  (2Wv2  +  20^)^/2 


tt)  +  a2  =  2, 


The  action  variables  defined  by 

J,  -  ^  Pjdq, 


give 


iL 


2no!i 


(4.88) 
(4.89) 


■2W 


(4.90) 

(4.91) 
(4.92) 
(4.93) 

(4.94) 
(4.95) 
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and 


2itZ2 


(^.96) 


Solving  Eqs.  (4.91),  (4-92),  (4.95)  and  (4.96)  for  a^    and  ^2  ^'^^   equating 
the  solutions 


and 


Hence 


"^^  JH  -  (Pj  -  2Hx^)/2  (^-97) 


2n 


■2W 


J  =  (p2  -  2Wv2) 


p2 
J  -  -  jt-sT^  (v^  +  -— i:L-) 
"*  (-2W) 


(4.98) 
(4.99) 
(4.100) 


The  form  of  these  latter  two  equations  suggests  that  two  new  variables 
be  defined 


a  =  |J.  + 


•2W 


M- 


(4.101) 


b  "  V  +  — Py   , 


so  that 


^ 


jtN/-2W  aa' 


Jy  »  -  nN/-2w  bb*  . 


(4.102) 

(4.103) 
(4.104) 


Both  a  and  b  are  very  similar  to  the  harmonic  oscillator  operators  in- 
troduced in  Chapter  2  and  by  analogy  one  would  expect  to  find  that 
quantities  proportional  to  ab*  and  a^b  are  also  constants  of  the  motion, 
This  is  in  fact  the  case.   Taking  real  and  imaginary  parts  of  these 
quantities  and  normalizing,  the  constants  are 
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K  - 


-^   (ab*  +  a*b) 


and 


L  =  -i  "LiW  (ab*  -  a*b). 
2 


(^.105) 


(4.106) 


The  third  independent  constant  can  be  chosen  to  be  the  difference  between 


i  and  J^ 


■2W 


(aa*  -  bb*) 


(4.107) 


It  Is  not  difficult  to  verify  that  these  constants  satisfy  the  commuta- 
tion relations  of  the  generators  of  SU2  or  R,  i.e., 


{K.  L)  =  -  2D 
(L,  D}  =  -  2K 
(D,  K)  =  -  2L  , 


where 


(a'^.  a}  -  (b*.  b)  = 


_  2i 


-2W 


(4.108) 
(4.109) 
(4.110) 

(4.111) 


It  should  be  noted  that  in  this  case  the  angle  variables  w^  and  w 
were  not  used.   This  occurred  only  because  a  similarity  was  noted  be- 
tween the  constants  of  the  motion  Jm  and  Jy  and  the  constants  of  the 
motion  for  the  harmonic  oscillator.   Implicit  in  the  definitions  of  a 
and  b  are  two  variables,  denoted  B  and  r,  which  are  conjugate  to  Jj^ 
Jy  but  these  are  not  the  angle  variables.   These  variables  are 


S  -  _  tan-1  .^_ 

In  si  -2\4x 


and 


I 
r  "  —  tan' 


■1   "^ 


(4.112) 


(4.113) 


2ff       ^/-2Wv 
t  is  easy  to  verify  that  these  are  conjugate  to  the  action  variables 
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for  (r,  Jy)  ■»  {S.  J^]   =  1  (4.1li*) 

and 

(r.  4x)  -  {s.  Jv)  •=  0  .  e+.ns) 

It  should  not  be  surprising  to  find  a  set  of  coordinates  conjugate  to 
the  action  variables  which  are  not  the  angle  variables  because  for  a 
given  set  of  momenta  the  set  of  canonical  coordinates  is  not  unique. 
This  occurs  because  one  can  always  add  a  total  time  derivative  to  the 
Lagrangian  which  is  a  function  of  q,  q  and  t. 

The  constants  K,  L,  and  D  have  remarl<ably  simple  interpretations. 
K  and  D  are  the  x  and  y  components  of  the  vector  P  =  R  •/-2W  where  R  is 
the  Runge  vector,  while  L  is  the  angular  momentum.   This  can  be  veri- 
fied by  substituting  Eqs.  (4.101)  and  (4.102)  into  Eqs.  (4.105)  through 
(4.107)  and  then  transforming  either  to  rectangular  or  polar  coordinates, 


CHAPTER  5 
SUMMARY 

First  of  all,  the  motivation  of  this  dissertation  was  to  find 
general  principles  under  which  symmetry  groups  of  Hamiltonlan  systems 
may  be  found.   One  is  guided  in  the  search  by  the  i<nown  facts  that  the 
Isotropic  harmonic  oscillator  and  the  Kepler  problems  have  such  symme- 
try groups,  as  has  been  shown,  for  example,  by  Saenz  [II],  who  reduced 
each  of  these  problems  to  force  free  motion  on  the  surface  of  a  hyper- 
sphere.   This  is  the  central  idea  behind  the  stereographic  parameters 
of  Laporte  and  Ralnich  [7];  Indeed  their  guiding  principle  is  that  for 
every  hypersphere  there  will  be  a  projection  associated  with  some  po- 
tential, which  will  have  accidental  degeneracy  on  that  account. 

However,  the  Ideas  expressed  In  this  dissertation  are  different 
and  depend  upon  two  stages  of  transformation.   The  first  Is  to  re- 
examine the  harmonic  oscillator,  using  another  mapping,  the  Hopf  map- 
ping, which  reduces  the  orbits  to  points.   Rotations  in  the  resulting 
space,  corresponding  to  SU_  transformations  of  the  original  space,  are 
readily  verified  to  form  a  symmetry  group.   The  exceptional  role  played 
by  the  harmonic  oscillator  in  this  regard  Is  simply  that  Its  Hamiltonlan 
Is  quadratic,  and  that  the  collection  of  all  quadratic  functions  over 
the  phase  space  forms  a  Lie  algebra  which  operates  on  phase  space  in 
a  particularly  Interesting  manner.   There  is  then  quite  generally  an 
Slip  group  commuting  with  any  normal  operator,  not  always  formed  from 
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quadratic  operators  except  in  the  isotropic  case,  for  which  the  normal 
operator  is  degenerate. 

The  second  stage  is  to  consider  what  Hamiitonians  are  equivalent 
to  quadratic  Hamiitonians  by  canonical  transformations.   In  this  re- 
gard one  can  say  almost  all  are.   Here,  one  reinterprets  the  Hamilton- 
Jacobi  equation.   If  it  is  thought  of  as  yielding  force-free  motion  in 
some  (generally  non-eucl idean)  phase  space,  then  a  further  transforma- 
tion to  polar  coordinates  will  yield  circular  motion  in  the  phase  planes, 
which  presumably  implies  a  quadratic  relation  between  the  coordinate  and 
conjugate  momentum. 

This  relation  is  delicate  on  two  accounts.   First  the  relationship 
in  one  phase  plane  may  depend  upon  separation  constants  arising  from 
other  phase  planes,  and  thus  irreconcilable  conflicts  may  arise  from 
Inverting  all  the  transformations  to  the  Hami I ton-Jacobi  system.   Secondly, 
with  an  eye  toward  quantum  mechanical  applications,  the  symmetry  operators 
may  well  be  transcendental  functions  of  the  coordinates  and  momenta,  mak- 
ing it  impossible  to  define  the  proper  quantum  mechanical  functions  of 
non-commuting  operators. 

The  organization  of  the  paper  was  as  follows:   In  the  first  chapter 
consideration  Is  given  to  the  Lie  algebra  generated  by  quadratic  functions 
over  the  phase  space  under  the  Poisson  bracket  operation,  and  its  relation 
to  different  homogeneous  function  spaces  also  taken  over  the  phase  space. 
It  is  seen  that  the  eigenfunct ions  of  these  quadratic  operators  in  the 
phase  space  itself  determine  the  structure  of  all  the  other  homogeneous 
spaces;  and  In  particular  there  is  a  simple  criterion  for  finding  a  set 
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of  operators  which  cormute  with  the  original  operator  under  Poisson 
brackets.   Namely,  it  is  that  they  be  constructed  from  products  of 
phase  space  eigenfunct ions  belonging  to  positive  and  negative  pairs 
of  eigenvalues.   These  correspond  to  the  "constants  of  the  motion" 
when  the  operator  Is  tal<en  to  be  the  Hamlltonlan.   It  is  also  seen 
that  such  positive  and  negative  pairs  will  exist  on  account  of  the 
general  existence  of  a  set  of  operators  formed  from  time-reversal  and 
coordinate  momentum  interchange  operators.   Even  though  the  constants 
will  form  a  linearly  independent  set  there  may  exist  functional  rela- 
tionships among  them  as  one  expects  to  find  at  most  (2n  -  1)  time  in- 
dependent constants  for  a  Hamlltonlan  with  n  degrees  of  freedom. 

The  constants  commuting  with  a  degenerate  normal  operator  generate 
an  SUn  group,  while  the  quadratic  functions  themselves  generate  an  Sp^, 
group  (symplectic  group  in  n  dimensions).   A  non-normal  operator  gen- 
erates a  somewhat  different  group  a  semidirect  product  [41]  of  a 

rotation  group  and  an  Abel i an  group,  while  normal,  non-degenerate 
operators  generate  an  SUp  group,  whose  operators  are  generally  not 
homogeneous  functions,  but  suitable  roots  of  the  latter. 

Secondly,  consideration  is  given  to  quadratic  Hamlltonlans  in  con- 
siderable detail,  the  second  chapter  being  devoted  to  isotropic  and 
anisotropic  harmonic  oscillators,  specifically  in  two  and  generally  in 
n  dimensions.   For  the  Isotropic  oscillators  a  particular  Hopf  mapping 
yields  a  direct  geometric  interpretation  of  the  constants  through  a 
canonical  mapping.   In  the  anisotropic  case  it  is  not  directly  a  geo- 
metrical transformation,  but  nevertheless  canonical  coordinates  and 
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momenta  are  found  such  that  the  momenta  are  n-1  of  the  generators  of 
the  SUn  group  which  generate  rotations  in  the  coordinate  planes. 

The  third  chapter  concerns  another  specific  class  of  quadratic 
Hami I tonians ,  which  may  be  interpreted  as  belonging  to  a  charged  par- 
ticle moving  in  a  uniform  magnetic  field  as  well  as  a  plane  harmonic 
oscillator  potential.   The  problem  is  of  some  interest  in  its  own  right, 
and  among  other  things,  it  is  learned  how  to  transfer  constants  of  the 
motion  to  a  rotating  coordinate  system.   In  particular  a  rule  is  given 
for  writing  Hamiltonians  with  central  potentials  when  a  uniform  magnetic 
field  Is  added  and  a  transformation  is  made  to  a  rotating  coordinate 
system  with  a  frequency  of  rotation  equal  to  the  Larmor  frequency.   The 
generators  of  the  symmetry  group  yield  a  set  of  intrinsic  coordinates, 
which  allow  a  clean  separation  of  the  magnetic  field  (or  angular  momen- 
tum) from  the  oscillatory  motion;  a  transformation  used  by  Goshen  and 
Lipl<in  [2k]    in  their  heuristic  treatment  of  rotational  bands.   In  addi- 
tion there  emerged  a  relation  between  the  oscillator  frequency,  the 
Larmor  frequency  and  a  parameter,  which  relation  yields  bounded  closed 
orbits  if  the  parameter  is  a  rational  number.   Finally,  by  taking  high 
field  and  low  field  limits,  one  obtains  Larmor  precession  of  the  har- 
monic oscillator  orbits  for  a  weak  field,  and  the  drift  of  the  cyclo- 
tron orbits  along  the  harmonic  oscillator  trajectories  in  the  high 
field  case.   Viewed  from  a  rotating  coordinate  system  with  the  Larmor 
frequency  all  orbits  become  pure  harmonic  oscillator  orbits.   Of  appli- 
cation to  general  problems  of  motion  in  a  magnetic  field  and  adiabatic 
constants,  this  problem  is  soluble  analytically,  so  that  the  magnetic 
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field  and  potential  field's  effects  may  be  seen  clearly.  A  careful 
treatment  is  given  of  the  known  constants  of  the  motion  of  the  cyclo- 
tron problem,  known  from  the  work  of  Johnson  and  Lippmann  [kl].   Par- 
ticular attention  is  paid  to  the  influence  of  the  vector  potential  on 
the  obvious  geometric  symmetry,  and  it  is  seen  that  canonical  rotations 
and  geometrical  rotations  yield  two  entirely  different  tyf>es  of  symme- 
try and  constants  of  the  motion. 

Although  the  specific  form  of  the  Hamiltonian  has  a  particularly 
rich  variety  of  interpretations  as  a  physical  problem  regarding  the 
motion  of  a  charged  particle,  mathematically  it  is  essentially  nothing 
but  a  (possibly  anisotropic)  harmonic  oscillator,  a  fact  which  reduces 
the  discovery  of  the  symmetry  group  to  a  previous  case,  and  renders 
the  analytic  solution  of  the  problem  trivial.   In  general  the  orbit  is 
a  Lissajous  figure  and  the  constants  of  the  motion  have  their  inter- 
pretation as  operators  which  move  the  boundaries  of  the  figure  or  change 
the  phase  of  the  orbit  in  various  ways,  preserving  the  sum  of  the  squares 
of  its  semiaxes  upon  which  the  energy  of  the  orbit  solely  depends. 

In  the  fourth  chapter  consideration  is  given  to  a  non-harmonic 
problem,  the  Kepler  problem,  which  is  well  known  to  be  accidentally  de- 
generate, and  equivalent  in  the  bound  orbits  at  least  to  the  harmonic 
oscillators.  As  a  new  result  it  is  found  that  there  is  an  SU3  symmetry 
group,  in  addition  to  the  known  R/^  symmetry  group,  the  additional  con- 
stants defining  the  line  of  nodes.   In  any  event  the  generators  of  the 
two  groups  are  functions  of  oneanother  and  not  subsets.   Thus,  one  has 
a  general  functional  relation  between  the  generators  of  these  two  groups. 
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As  a  general  rule  one  can  say  that  any  Hamiltonian  which  is  separable 
will  possess  an   SU^  symmetry  even  though  the  frequencies  in  each  angle 
coordinate  may  not  be  equal. 

Since  it  is  of  considerable  importance  for  the  further  application 
of  this  technique,  the  relation  of  the  SUo  group  arising  from  different 
separations  of  the  Ham i 1  ton -Jacob i  equation  in  polar  and  parabolic  co- 
ordinates is  studied.   It  is  found  that  considerable  attention  has  to 
be  paid  to  the  problem  of  finding  rational  constants  of  the  motion. 

This  work  offers  many  prospects  for  continuation.   There  are   still 
many  potentials,  for  which  hope  of  finding  rational  constants  seems 
slim,  but  which  nevertheless  ought  to  be  investigated.   Among  these  are 
the  Stark  effect  (Redmond  recently  found  a  new  constant  [kS]) ,    the 
problem  of  two  fixed  attracting  centers  [46],  the  motion  of  a  charged 
particle  in  the  field  of  a  magnetic  dipole,  as  well  as  the  relativistic 
Kepler  problem  which  Biedenharn  has  investigated  recently  [27,  28].   In 
addition,  the  relations  uncovered  between  the  orthogonal  groups  and 
unitary  groups  might  lead  to  some  interesting  formulas  for  the  higher 
dimensional  groups,  but  a  point  of  particular  interest  is  the  appear- 
ance of  the  angular  momentum  operators  as  a  subgroup  of  SU,.   Since 
they  depend  on  two  of  the  three  coordinates  of  phase  space,  one  can 
exploit  the  symmetry  of  SU^  as  a  permutation  of  its  generators  to  ob- 
tain SU2  groups  involving  the  remaining  coordinate,  and  hence  the  ra- 
dial action.   In  quantum  mechanical  terms,  this  yields  another  view- 
point toward  the  known  fact  that  coulomb  integrals  behave  as  Wlgner 
coefficients  [47]. 
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In  conclusion  it  should  be  pointed  out  that  due  to  the  nonlinear- 
Ity  of  the  various  relations  found,  their  extension  to  quantum  mechani 
cal  operators  will  require  exceptionally  careful  scrutiny. 
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